-v* »»«*»•;«*=? ^ «i 


AD-A0C8  847 

IONOSPHERIC  MODIFICATION:  ITS  EFFETS 

ON  RADAR  AND  SATELLITE  COMMUNICATION 
SYSTEMS 

G.  F.  Var.Blaricum  , et  al 
General  Research  Corporation 


Prepare  d for : 

Army  Missile  Command 

Defense  Advanced  Research  Projects  Agency 


October  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


. * \ -^»3S^^^f®V‘/i*iSOT^I 


ARP  A Order  Number: 

Program  Code  Number: 

Short  Title  of  Work: 

Name  of  Contractor: 

Amount  of  Contract: 

Contract  Number: 

Effective  Date  of  Contract: 
Contract  Expiration  Date: 
Project  Engineer: 

Phone  Number: 


958,  Amendment  No.  5 
4E20 

STO  Supporting  Studies: 

Signal  Transmission  Data  Analysis 
General  Research  Corporation 
$36,829  (Task  18) 

DAAH01-7  3-C-0310 
31  May  1974 
31  December  1974 
A.  D.  Stathacopoulos 
(805)  964-7724,  Ext.  412 


j Met'  A V, 

I ’ 1 


f; 


In  addition  to  approval  by  the  Project  Leader 
and  Department  Head,  General  Research  Corporation 
reports  ar$  subject  to  independent  review  by 
'staff" tdfemb^  not  connected  with  the  project. 

- This  repoj^was  reviewed  by  P.  J.  Redmond. 


ThtFYiews  and  conclusions  contained  in  this  document  are  those  of  the 
humlrs  and  should  not  be  interpreted  as  necessarily  representing  the 
pffftrai  policies,  either  expressed  or  implied,  of  the  Defense  Advanced 
ke§*OTCh' PrdjertS'TSgincy  or  the  U.S.  Government, 


This  research  was  supported  by  the  Defense  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by  the  U.S.  Army 
Missile  Command  under  Contract  DAAH01-73-C-0310. 


L 


§ 

I 


security  classification  of  this  page  jhtimi  lute  F.nterv.h 


1.  REPORT  NUMBER 

jRM-1911 


REPORT  DOCUMENTATION  PAGE 


2.  60VT  ACCESSION  NO 


4.  TITLE  (mnd  Submit) 

Ionospheric  Modification:  Its  Effects  on  Radar  and 

Satellite  Communication  Systems 


7 AUTHORS; 

G.  F.  VanBlaricum,  H.  S.  Ostrowsky 


*■  PERFORMING  ORGANIZATION  name  and  address 

General  Research  Corporation 
P.0.  Box  3587 
Santa  Barbara,  CA  93105 


II.  CONTROLLING  office  name  and  address 
Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  VA  22209 


M monitoring  AGENCY  NAME  4 ADDRESSfi/  dtllerent  Item  Controlling  Other) 

US  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 
Attn:  AMSMI-RNS 


HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 RECIPIENT’S  CATALOG  NUMBER 
S.  TYPE  OF  REPORT  a PERIOD  COVEI 


COVERED 

Releasable  Memorandum 
31  May  - 30  September  1974 


(.  PERFORMING  ORG.  REPORT  NUMBER 

RM-1911 


1.  CONTRACT  OR  GRANT  kiumBERC') 

DAAH01-7  3-C-0310 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

ARPA  Order  No.  958 
Amendment  No.  5 


12-  REPORT  DATE 

October  1974 


13  NUMBER  OF  PAGCS 


IS.  SECURITY  CLASS,  (ol  (his  report) 

UNCLASSIFIED 


IS«.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (ot  the  ebatrect  entered  in  Block  20.  it  different  trom  Bepar t) 


ib  supplementary  notes 


IS  KEY  WOROS  ( Cocir tnue  on  reveree  tide  it  neceasmry  end  identity  by  b/ocA  number) 

Scini'.illation  . . 

Reproduced  by 

Ionospheric  modification  NATIONAL  TECHNICAL 

Radar  systems  INFORMATION  SERVICE 

Satellite  communication  systems  USSFHnXw"'vA.Ci2isi,c* 


cO  ABSTRACT  (C-mfinu#  on  revtt%*  tide  1/  necMiery  mnd  identity  bv  block  number < 

High-power  RF  radio  transmissions  have  been  used  to  create  significant  local  dis- 
turbances of  the  temperature  and  electron  density  in  the  F-region  of  the 
ionosphere.  This  artificial  modification  of  the  ionosphere  can  create  appreciable 
levels  of  amplitude  and  phase  scintillation  in  radio  signals  passing  through  the 
disturbed  region.  A number  of  measurements  of  such  scintillation  were  conducted 
by  Aeronomy  Corporation  in  an  experimental  program  known  as  PRAIRIE  SMOKE. 

(Continued) 


DD  1 j an  73  1473  EDITION  OF  I NOV  AS  IS  OBSOLETE  • 

I 


• security  classification  of  this  page  (iwien  f> atm  Fm ere/> 


SECURITY  CLASSIFICATION  or  THIS  PAOEfWhtw  Umf  Knitted)  

(Block  20  continued) 

This  report  describes  and  reorganizes  some  of  the  PRAIRIE  SMOKE  experimental 
data  and  presents  a brief  overview  of  scintillation  theory.  The  experimental 
data  were  used  to  develop  a model  for  predicting  scintillation  effects  along 
any  particular  line-of-sight  through  the  disturbed  region.  That  model  is 
presented  here  and  predictions  made  using  it  are  compared  to  some  of  the  experi- 
mental data. 

Scintillation  can  degrade  radar  detection  probability,  radar  tracking  accuracy, 
and  satellite-to-ground  communication  error  rates.  These  degradation  were 
computed  parametrically  as  functions  of  the  scintillation  index.  The  results 
presented  in  this  report  can  be  used  to  determine  the  degradation) produced  by 
natural  or  nuclear-induced  scintillations  as  well  as  scintillations  produced 
by  RF  heating  of  the  ionosphere. 


\<L, 


SECURITY  Cl  ASSiriCATlON  OF  THIS  PAGE'WhYn  D.1  in  Fnirmti 


SI  ;'J?.',-t5^V ?£?•**<.  J ^,1,^  !<T^  jII’VAj  ' -,  MIT  t!':'  Hi  '■• . j 'rr  “''1^?^^^.?.1’.V^/^T7VTr»^r-;^r,;j|rjp»l~_  svJr- 


ACKNOWLEDGMENT 

The  authors  thank  Dr.  S.  A.  Bowhill,  Dr.  E.  K.  Walton,  and 
Mr.  D.  R.  Ward  of  Aeronomy  Corporation  for  their  cooperation  in  providing 
the  data  from  Prairie  Smoke,  for  their  helpful  discussions  regarding  the 
experimental  program,  and  for  their  suggestions  of  the  appropriate  form 
for  the  scintillation  model. 


sr.*;s  rei.-i-.fe 


I 


K: 

*\ 

jfel 


ABSTRACT 


a'3-  ? ^Tf'T-iWJ  7T<?W  WT^rr^Mflr-.T-Tc.-i " 


High-power  RF  radio  transmissions  have  been  used  to  create  significant 
local  disturbances  of  the  temperature  and  electron  density  in  the  F-region 
of  the  ionosphere.  This  artificial  modification  of  the  ionosphere  can 
create  appreciable  levels  of  amplitude  and  phase  scintillation  in  radio 
signals  passing  through  the  disturbed  region.  A number  of  measurements 
of  such  scintillation  were  conducted  by  Aeronomy  Corporation  in  an 
experimental  program  known  as  PRAIRIE  SMOKE. 


This  report  describes  and  reorganizes  some  of  the  PRAIRIE  SMOKE 
experimental  data  and  presents  a brief  overview  of  scintillation  theory. 

The  experimental  data  were  used  to  develop  a model  for  predicting  scintilla- 
tion effects  along  any  particular  line-of-sight  through  the  disturbed  re- 
gion. That  model  is  presented  here  and  predictions  made  using  it  are  com- 
pared to  some  of  the  experimental  data. 


Scintillation  can  degrade  radar  detection  probability,  radar 
tracking  accuracy,  and  satellite-to-ground  communication  error  rates. 

These  degradations  were  computed  parametrically  as  functions  of  the 
scintillation  index.  The  results  presented  in  this  report  can  be  used  to 
determine  the  degradation  produced  by  natural  or  nuclear- induced  scintilla- 
tions as  well  as  scintillations  produced  by  RF  heating  of  the  ionosphere. 
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INTRODUCTION  AND  SUMMARY 


Since  about  1970  experimenters  have  been  heating  the  F-region  of 
the  ionosphere  with  high-power , ground-based  RF  transmitters,  creating 
local  perturbations  of  up  to  30  percent  in  the  temperature  and  1 percent 
in  the  electron  density  of  the  ionospheric  plasma.  These  artificially 
produced  disturbances  are  observable  in  modified  optical  emissions,  in 
radar  backscatter  and  forward-scatter,  in  lonog  .ms,  and  in  scintillations 
produced  in  radio  signals  traversing  the  disturbed  region.^  Although 
artificial  modification  of  the  ionosphere  has  been  carried  out  in  the 
past  by  release  of  chemicals  from  rockets  and  detonation  of  atomic  bombs, 
RF  heating  of  the  ionosphere  is,  in  contrast  to  these  means,  controllable, 
repeatable,  and  rapidly  self-reversing. 

The  ionospheric  modification  experiments  have  been  conducted  using 
a transmitter  constructed  near  Platteville,  Colorado,  by  the  Institute 
of  Telecommunication  Sciences  of  the  US  Department  of  Commerce.  (A  few 
experiments  have  also  been  run  using  the  300-meter  antenna  at  Arecibo  Ob- 
servatory, Puerto  Rico.)  Project  IVORY  CORAL,  sponsored  by  the  Defense 
Advanced  Research  Projects  Agency  (DARPA) , has  encompassed  a number  of 
theoretical  and  laboratory  studies  conducted  to  develop  an  understanding 
of  ionospheric  modification,  and  a number  of  diagnostic  field  experiments 
(known  as  Project  PRAIRIE  SMOKE)  were  conducted  with  the  Platteville 
transmitter  during  1972-73.  Studies  of  possible  systems  applications  of 
ionospheric  modification — for  example,  as  a communications  medium — and 
some  feasibility  demonstrations  were  also  sponsored  by  DARPA  as  Project 
SMOOTH  LEDGE. 

The  present  study,  a part  of  Project  SMOOTH  LEDGE,  was  undertaken 
to  examine  the  possible  systems  implications  of  the  scintillation 
produced  in  radio  signals  propagating  through  an  artificially  modified 
ionosphere.  The  objectives  of  this  study  were  to: 

• Describe  and  organize  the  scintillation  data  measured  in  the 
PRAIRIE  SMOKE  experiments 
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• Fit  the  measured  data  to  a theoretical  scintillation  model 
which  car.  then  be  used  to  predict  further  effects 

• Relate  the  scintillation  data  to  parameters  which  describe 

or  determine  the  degradation  of  radar  or  communication  system 
performance 

As  part  of  the  PRAIRIE  SMOKE  experiments,  Aeronomy  Corporation 
measured  the  scintillations  produced  by  the  modified  ionosphere  in  VHF 
signals  transmitted  from  polar  orbiting  and  geostationary  satellites  to 
ground  receiving  stations.  Aeronomy  Corporation  used  their  data  to 
diagnose  the  size  and  shape  of  the  modified  region  and  the  ionization 
levels  and  correlation  lengths  which  determine  its  small-scale  structure. 
Much  of  their  measured  data,  along  with  descriptions  of  their  experimental 
configurations  and  the  satellite  systems  employed,  were  made  available  to 
us. 


The  measured  scintillations  due  to  the  artificially  modified 

ionosphere  were  compared  to  the  predictions  of  the  Briggs  and  Parkin  model 

2 

for  scintillation  in  the  natural  ionosphere.  Although  a Briggs  and  Parkin 
model  had  been  assumed  in  deriving  many  of  the  parameters  describing  the 
modified  region,  that  modal  failed  to  predict  the  observed  variation  in 
scintillation  with  the  angle  between  the  geomagnetic  field  and  the  trans- 
mitter-to-receiver  line-of-sight.  The  disagreement  was  soon  attributed 
to  the  finite  dimensions  of  the  modified  region  as  compared  to  the  assumed 
infinite  horizontal  extent  in  the  Priggs  and  Parkin  model.  When  the  model 
was  refined  to  account  for  the  finite  region,  good  agreement  between 
measurements  and  model  predictions  was  obtained. 

The  principal  effects  of  the  artificially  modified  ionosphere  on  a 
radar  forced  to  look  through  it  are  amplitude  scintillations  in  the 

ft 

received  signal  and  scintillations  in  angle  of  arrival  of  the  signal. 


Phase  coherence  in  the  signal  may  be  degraded  as  well,  but  the  effects 
should  be  negligibly  small  over  the  short  duration  of  most  radar  pulses. 
Effects  of  phase  coherence  have  not  been  treated  explicitly  in  this  study. 
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For  an  assumed  Rice  probability  distribution  of  the  amplitude  fluctuations, 
the  degradation  In  radar  detection  probability  and  in  detection  range  were 
determined  parametrically  for  a wide  range  of  slgnal-to-noise  ratios  and 
scintillation  indexes.  The  radar  angle  errors  produced  by  various  levels 
of  scintillation  were  also  developed.  In  order  to  assess  the  overall 
effect  of  the  artifically  modified  ionosphere  on  radar  tracking,  a Monte 
Carlo  computer  simulation  was  used  to  determine  the  resultant  errors  in 
predicting  the  intercept  point  for  a ballistic  missile  in  two'  cases: 
first,  when  the  radar  viewed  the  missile  through  the  modified  ionosphere 
just  as  the  missile  appeared  above  the  horizon,  and,  second,  when  the 
disturbed  tracking  occurred  only  at  the  end  of  the  tracking  interval. 

The  resultant  intercept  point  errors  were  computed  for  many  choices  of 
scintillation  index  and  parameters  of  the  modified  region. 

The  degradation  in  the  bit  error  rate  in  satellite-to-ground  DPSK 
(differential  phase-shift  keying)  and  NFSK  (non-coherent  frequency-shift 
keying)  digital  communications  was  computed  for  the  assumed  Rice  distri- 
bution of  amplitude  fluctuations  as  a function  of  the  scintillation  index. 

Most  of  the  parametric  results  derived  for  radar  detection  and 
tracking  and  for  satellite  communication  depend  only  upon  the  intensity 
of  the  scintillation  and  not  its  source,  sc  that  the  results  could  be 
applied  to  natural  scintillation  or  to  nuclear-weapon-induced  scintilla- 
tions as  well  as  to  those  produced  by  ionospheric  heating. 

Briefly  summarized,  the  principal  conclusions  of  the  study  are: 

• RF  heating  of  the  ionosphere  can  produce  scintillations  about 
twice  as  great  as  those  observed  in  the  normal  daytime  iono- 
sphere in  the  middle  latitudes.  The  largest  amplitude  scin- 
tillation measured  in  PRAIRIE  SMOKE  at  150  MHz  was  about  10 
dB  peak-to-peak;  most  values  were  much  less. 


• The  enhanced  scintillation  is  produced  only  along  lines  of 
sight  through  the  heated  region,  which,  for  the  Platteville 
heater,  subtends  at  most  about  a 20-degree  cone  angle  at 
ground  receiving  sites.  The  largest  scintillations  occur 
when  the  line-of-sight  is  parallel  to  the  geomagnetic  field. 

• Amplitude  and  phase  scintillations  can  degrade  radar  detec- 
tion and  tracking  and  satellite  communication  error  rates. 
However,  possible  intentional  degradations  produced  by  iono- 
spheric modification  can  be  largely  discounted  because: 

1.  The  scintillations  are  produced  only  in  a small  angular 

sector  of  the  sky. 

2.  Scintillation  scales  with  frequency  at  a rate  between 

-1  -2 

f and  f . Systems  operating  at  UHF  and  above 
should  experience  only  a small  fraction  of  the  scintil- 
lation measured  in  the  experiments  at  150  MHz. 

3.  The  margin  usually  engineered  into  systems  to  account 
for  severe  natural  scintillation  may  offset  the  worst 
effects  of  artificially  produced  scintillation.  At  any 
rate,  the  largest  natural  scintillations  are  likely  to 
be  several  times  worse  than  those  produced  by  RF  heat- 
ing. (For  example,  8 dB  fading  has  been  observed  on 
satellite  communication  channels  at  6 GHz.) 

Section  2 of  this  report  describes  the  natural  Ionosphere  and  iono- 
spheric modification.  Section  3 summarizes  a simple  theory  of  ionospheric 
scintillation  and  describes  the  extent  of  natural  scintillation.  Section 
4 describes  the  satellite  transmission  experiments  in  PRAIRIE  SMOKE  and 
compares  the  measured  data  with  values  derived  from  the  simple  scintilla- 
tion model;  Sec.  5 presents  the  refined  scintillation  model. 

The  degradation  of  radar  detection  and  tracking  is  discussed  in  Sec. 
6 and  the  degradation  of  satellite  communication  is  discussed  in  Sec.  7. 
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2 IONOSPHERIC  MODIFICATION 

2.1  THE  NATURAL  IONOSPHERE 

The  Ionosphere  is  the  portion  of  the  earth's  upper  atmosphere  in 
which  enough  of  the  constituent  gases  are  ionized  to  affect  the  propaga- 
tion of  radio  waves.  For  practical  purposes  the  ionosphere  extends  from 
about  50  km  altitude  out  to  several  earth  radii.  It  is  conventionally 
divided  into  three  regions:  D (50  to  90  km),  E (90  to  130  km),  and  F 

(above  130  km).  In  addition,  the  E-  and  F-regions  are  divided  into 
layers,  of  which  the  F2  layer,  at  about  300  km,  is  the  most  important 
from  the  point-of-view  of  radio  propagation. 

At  300  km  altitude,  1 cubic  meter  of  the  ionosphere  contains  about 
14  15 

10  to  10  neutral  particles  (CIRA  1965  model  atmosphere)  and  about 
10^*  to  10"^  electrons  and  an  equal  number  of  positive  ions,  chiefly  0+. 
Because  of  their  much  smaller  masses,  the  electrons  are  more  mobile  than 
the  positive  ions  and  hence  determine  most  of  the  electrical  properties 
of  the  ionosphere.  The  electrons  oscillate  weakly  at  the  plasma  frequency, 
which  is  proportional  to  the  square  root  of  the  electron  density,  and 
reflect  or  absorb  radio  waves  at  that  frequency  or  lower.  The  highest 
plasma  frequency  in  the  ionosphere  is  called  the  critical  frequency  of 
the  F-layer;  it  typically  ranges  from  !*  to  10  MHz. 

Because  of  the  earth's  magnetic  field,  the  ionosphere  is  birefrin- 
gent  with  two  modes  of  propagation  which  travel  at  different  velocities 
and  which  refract  differently.  The  modes  are  usually  denoted  "0,"  for 
ordinary,  and  "X,"  for  extraordinary.  The  0-mode  has  approximately 
right-circular  polarization  and  tends  to  deviate  normal  to  the  magnetic 
field;  the  X-mode  is  approximately  left-circular  and  tends  to  deviate 
parallel  to  the  magnetic  field.  In  practice,  the  0-mode  and  the  X-mode 
have  different  critical  frequencies,  denoted  fQ  and  f^  . Since  there 
are  critical  frequencies  associated  with  all  of  the  ionospheric  layers, 
these  frequencies  in  the  F2  layer  are  usually  labeled  fQF2  and  fxF2. 


5 


2.2  MODIFYING  THE  ATMOSPHERE 

Local  perturbations  of  the  Ionosphere  ar,i  being  created  by  deliver- 
ing large  amounts  of  energy  into  the  F2  layer  from  a high  power  RF  trans- 
mitter operating  near  the  critical  frequency.  The  transmitter,  located 
at  Platteville,  Colorado  (40.18®N,  104.73°W),  has  a peak  power  of  about 
2 MW.  The  antenna,  designed  to  operate  over  5 to  10  MHz,  is  a ring  array 
with  about  a 15-degree  half-power  beunwidth  that  can  be  driven  to  generate 
right-circular  (<vO-mode)  and  left-circular  (^X-mode)  polarized  waves. 

The  transmitter  and  antenna  yield  a power-aperture  product  of  about 
10  2 

10  W#m  (the  300-meter  antenna  at  Arecibo  Observatory,  Puerto  Rico, 
has  also  been  employed  in  ionospheric  modification  experiments).^ 

A detailed  discussion  of  the  physical  processes  involved  in  iono- 
spheric heating  is  beyond  the  scope  of  this  paper  (for  such  a discussion, 
the  reader  is  referred  to  Ref.  3).  However,  a very  simplified  descrip- 
tion is  as  follows.:  ■ , 


Near  the  altitude  where  the  plasma  frequency  approximates  the  trans- 
mitted frequency,  the  electrons  are  driven  into  oscillation  and  convert 
some  of  the  RF  energy  into  random  thermal  motion  as  they  collide  with 
ions  and  neutral  particles.  Because  of  this  increased  thermal  motion, 
the  ambient  electron  temperature  is  increased.  Heat  is  conducted  away 
from  the  altitude  of  absorption  by  electrons  diffusing  along  the  geomag- 
netic field  lines.  Because  of  instabilities  the  diffusion  process  results 
in  a striated  structure  of  temperature  irregularities  extended  along  the 
field  lines.  This  diffusion  also  modifies  the  electron  density  in  the 
heated  region,  resulting  in  local  variations  in  the  mean  electron  density. 
The  local  variation  in  the  electron  density  is  important  in  determining 
the  diffraction  and  backscattering  properties  of  the  heated  region.  A 
schematic  view  of  ionospheric  heating  is  shown  in  Fig.  2.1. 

2.3  EXPERIMENTAL  STUDIES  OF  THE  MODIFIED  IONOSPHERE 

Most  of  the  experimental  techniques  used  to  study  the  normal  iono- 
sphere can  be  employed  to  measure  the  properties  of  the  modified 
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Figure  2.1.  Schematic  of  Ionospheric  Modification 


ionosphere.  One  of  the  principal  tools  used  in  ionospheric  studies  is 
the  ionosonde,  a kind  of  swept-f requency  radar  used  to  produce  plots — 
called  ionograms — of  the  height  of  radio  reflection  as  a function  of  fre- 
quency and  to  determine,  among  other  properties,  the  critical  frequencies 
of  each  layer.  Ionograms  taken  of  the  normal  ionosphere  usually  show 
well-defined  narrow  traces  of  reflection  height  in  the  F-region,  indicat- 
ing that  reflections  occur  at  one  height  for  each  frequency.  Ionograms 
of  the  heated  region,  however,  show  the  occurrence  of  "spread-F,"  a widen- 
ing of 'the  trace  indicating  scattering  from  irregularities  over  a con- 

4 

siderable  thickness  of  the  ionosphere.  Spread-F  occurs  naturally,  but 
usually  only  at  night  (at  Boulder),  and  is  usually  accompanied  by  increased 
scintillation.  The  heating  of  the  ionosphere  creates  spread-F  irregulari- 
ties even  in  daytime  shortly  after  the  transmitter  is  turned  on. 

Modifications  in  the  ionosphere  have  been  observed  in  optical  emis- 
7 8 

sions.  * Emissions  at  630  nm  are  normally  produced  in  the  dissociative 
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recombination  of  electrons  and  0.+  ions,  a reaction  which  proceeds  more 

1 7 

slowly  as  the  temperature  is  raised.  Biondi,  et  al.,  measured  a decrease 
in  630-nm  emissions  in  the  heated  ionosphere  corresponding  to  about  a 35- 

g 

percent  increase  in  the  temperature.  Evans » et  al.,  reported  an  enhance- 
ment in  the  1.27-ym  emission  of  molecular  oxygen  believed  also  to  be  due 
to  the  increased  temperature  in  the  heated  region. 


| Measurements  of  radar  backseat ter  were  used  to  determine  the  radar 

cross  section  of  the  heated  region,  the  frequency  dependence  of  that  cross 
section,  and  to  investigate  the  aspect  sensitivity  of  the  backscatter. 

If  the  irregularities  (or  striations)  in  the  heated  region  are  highly 
elongated  and  aligned  with  the  geomagnetic  field,  reflections  would  be 
expected  to  be  very  sensitive  to  aspect  angle.  A radio  wave  Incident  on 
the  Irregularities  making  a small  angle  6 to  the  normal  plane  to  the 
field  line  would  be  scattered  in  a cone  at  angle  -0  to  the  field  normal. 

f ■ 

F-  A simple  analogy  is  the  scattering  of  a flashlight  beam  by  a bundle  of 

parallel  glass  rods.  The  predicted  aspect  sensitivity — confirming  the 
field  alignment  of  the  Irregularities — was  observed  by  bistatic  measure- 
ments from  sites  about  350  km  (magnetic)  south  of  Platteville  and  by  RAM 
radar  measurements  from  White  Sands  Missile  Range  somewhat  further  south. 

A fourth  type  of  experiment  measured  the  scintillation  statistics 
of  satellite  radio  signals  passing  through  the  heated  region  to  receiving 
stations  on  the  ground.  The  experiments  were  carried  out  by  Aeronomy 
Corporation  during  PRAIRIE  SMOKE  I-V  using  both  geostationary  and  polar- 
orbiting  satellites  transmitting  in  the  VHF  and  UHF  bands.  Ground  receiv- 
ing sites  were  chosen  to  view  the  satellites  along  lines-of-sight  passing 
through  the  heated  region.  In  viewing  the  geostationary  satellite,  the 
llne-of-sight  remained  fixed,  but  in  viewing  the  orbiting  satellite  the 
line-of-sight  swept  through  the  heated  region  in  a few  minutes  as  the 
satellite  made  a north-to-south  (or  south-to-north)  orbital  pass.  For 
most  of  the  orbital  experiments  the  receiving  stations  were  placed  near 
the  down-field  point  (about  120  km  (magnetic)  north  of  Platteville)  so 


r^i(i«,ri 


that  the  line-of-sight  through  the  heated  region  would  be  as  nearly  paral 
lei  to  the  magnetic  field  as  possible. 

The  satellite  transmission  experiments  provided  considerable  infor- 
mation on  the  size  and  structure  of  the  heated  region  and  Its  scattering 
properties,  including: 

• Field  alignment  of  irregularities 

• Dimensions  and  orientation  of  correlation  ellipsoids  for 
large  and  small  irregularities 

• Drift  velocity  of  irregularities 

• Altitude  and  horizontal  dependence  of  signal  fluctuations 

• Intensity,  time  scale,  and  frequency  dependence  of  signal 
fluctuations 

Since  ionospheric  scintillation  induced  by  ionospheric  modification 
is  the  principal  focus  of  this  paper,  the  results  of  the  satellite  experi 
ments  are  very  important  tc  our  conclusions.  Most  of  the  results  of  the 
satellite  transmission  experiments  require  interpretation  in  terms  of  the 
theory  of  ionospheric  scintillation.  In  order  to  make  those  interpreta- 
tions most  readily  understandable,  we  defer  further  discussion  of  the 
satellite  experiments  to  Sec.  4 and  present  some  of  the  most  important 
aspects  of  scintillation  theory  in  Sec.  3. 
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THEORY  OF  IONOSPHERIC  SCINTILLATION 


Shortly  after  the  existence  of  radio  stars  and  galaxies  was  dis- 
covered, fluctuations  In  the  amplitude  of  the  received  signals  were  ob- 
served. By  studying  the  correlation  between  fluctuations  observed  at 
separated  antennas  It  was  determined  that  the  signal  scintillations  are 
primarily  produced  in  the  ionosphere,  and  principally  in  the  F-region 
where  the  ionization  is  strongest.  The  advent  of  artificial  earth  satel- 
lites in  1957  led  to  increased  interest  in  ionospheric  scintillation,  both 
to  understand  possible  degradation  of  communication  signals  as  well  as  to 
employ  satellites  as  tools  in  developing  an  understanding  of  che  structure 
and  propagation  properties  of  the  ionosphere. 

Scintillations  induced  in  radio  waves  passing  through  the  ionosphere 
have  been  studied  theoretically  for  about  twenty-five  years.  *'  Most 
of  the  results  of  interest  to  the  present  problem  are  drawn  from  Refs. 

2,  10,  and  11. 

3.1  THE  IONOSPHERE  AS  A PEASE  SCREEN 

A common  approach  in  these  papers  is  to  equate  the  effect  of  the 

ionosphere  to  that  of  a thin  screen  which  modulates  the  phase  or  amplitude 

of  the  wave  passing  through  it.  (Bowhill^  worked  with  an  amplitude- 

modulating  screen,  but  had  shown  earlier  the  equivalence  between  shallow 

amplitude  and  phase  screens.  We  shall  consider  only  phase  screens  here.) 

The  statistics  of  the  phase  screen  are  described  by  a two-dimensional 

stationary  covariance  function.  Bowhill  considered  only  shallow  screens 

whose  RMS  phase  variation  was  less  than  1 radian;  Mercier^  considered 

2 

large  phase  variations;  and  Briggs  and  Parkin  extended  some  of  Mercier's 
results  to  the  non-isotropic  phase  screen. 

The  complex  amplitude  of  a signal  propagating  in  the  positive  z- 
dlrection  emerging  from  the  random  phase  screen  at  the  plane  z ■ 0 (see 
Fig.  3.1)  is 


A(x,y,0)  - Aq  exp (j4> (x,y)j 
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WAVE 


Figure  3.1.  Diffraction  by  a Random  Phase  Screen 


where  $(x,y)  Is  * real  stationary  random  function  with  covariance 
^(Xi  " *j»  y±  " yj)  " E[<Kx1,y1)<j>(Xj,yJ)]  » where  E[ • ] denotes  ex- 
pected value  or  ensemble  average.  The  incident  wave  is  assumed  to  be 

plane  at  the  phase  screen.  An  observer  in  a plane  removed  from  the  phase 

1/2 

screen  observes  the  amplitude  R ■ (AA*)  of  the  diffracted  and  re- 
fracted wave,  where  A(x,y,z)  is  the  complex  amplitude.  Our  concern  is 
with  the  statistics  of  R and  its  covariance  function  in  the  observation 
plane. 

3.1.1  The  Phase  Screen — Briggs  and  Parkin  Theory 

Because  of  the  ionization  in  the  F-region  of  the  ionosphere,  the 
refractive  index  deviates  from  unity  and  the  total  phase  change  of  a wave 
propagating  through  the  region  differs  from  the  phase  change  accumulated 
over  an  equal  path  in  free  space.  At  frequencies  well  above  the  plasma 
frequency  of  the  medium,  typically  of  the  order  of  5 MHz,  the  refractive 
index  is  given  approximately  by 
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V “ 1 - XV  N/(2ir) 

e 


where  X is  the  wavelength,  r is  the  classical  electron  radius 


/ 2 2 —15  i 

(■  e /(Aireonic  ) “ 2.82  x 10  m)  , and  N is  the  electron  density. 


The  deviation  of  the  refractive  index  from  the  mean  value  due  to  local 


variations  in  the  electron  density  is  given  by  Ay  ■ -X  reAN/(2ir)  , where 


AN  is  the  deviation  of  N from  its  mean  value.  Then  the  phase  change, 
in  traveling  a distance  L in  the  medium,  relative  to  the  nhhse  change 
In  traveling  an  equal  distance  in  a uniform  medium  with  the  same  mean 
electron  density,  is 


2tt 

X 


Ld 

I 


Ay  dJt 


-r  X 
e 


i-i 
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AN  d£ 


Briggs  and  Parkin  derived  expressions  for  the  RMS  phase  fluctua- 
tions in  a wave  passing  through  a slab  in  which  the  electron  density  varia 
tlon  consists  of  a large  number  of  strlations  elongated  along  the  direc- 
tion of  the  geomagnetic  field.  The  correlation  of  the  electron  density 
in  the  medium  was  taken  to  be 


PN(r,s)  - exp|-r2/r2  - s^aro)2] 


where  r and  s are  cylindrical  coordinates  such  that  s *s  measured 
along  the  magnetic  field  and  r is  measured  perpendicular  to  it  (see 
Fig.  3.2);  tq  is  a characteristic  dimension  of  the  medium  which  is  very 
important  in  the  ensuing  discussion. 

For  a wave  incident  at  angle  i upon  a slab  of  thickness  Ah  , the 
path  length  through  the  slab  is  L » Ah  sec  i (see  Fig.  3.3).  Briggs 
and  tarkin  derived  the  mean  square  phase  variation  for  the  wave  to  be 
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Figure  3.2.  Dimensions  of  Field-Aligned  Striations 
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E(*2) 


ff^2r2X2E(AN2)Ahotro  sec  1 
(a2  sin2  ip  + cos2  ip)1/2 


(3.1) 


where  is  the  angle  between  the  magnetic  field  vector  B and  the  nor- 
mal to  the  incident  wave.  The  correlation  function  for  the  phase  fluctua- 
tion in  the  wave  emerging  from  the  slab  is 


P^u.v) 


„ 2,2  2,,0  .21 
exp  -u  /r  - v / (0r  ) 
o o 


(3.2) 


where 


6 


sin2ij> 


+ cos  2(|<)^2 


(3.3) 
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3.2  STATISTICS  OF  THE  DIFFRACTED  WAVE 

As  tha  wava  propagates  away  from  the  phase  screen,  the  diffracted 
components  of  the  wave  begin  to  Interfere  constructively  and  destructively 
with  the  undiffracted  or  "specular"  components  so  that  the  amplitude  of 
the  signal  scintillates.  In  this  section  we  are  concerned  with  the  fol- 
lowing Issues: 

• What  are  the  statistical  properties  of  the  amplitude  and 

phase  scintillations  of  the  observed  wave? 

\ 

• How  do  the  statistics  vary  with  the  distance  z between 

o 

the  phase  plane  and  the  observer  plane? 

• How  are  the  scintillations  correlated  In  the  observer  plane 
and  how  does  the  correlation  relate  to  the  correlation  In 
the  phase  plan?? 

The  measurements  In  the  observer  plane  are  usually  either  of 
1/2  2 

R ■ (AA*)  or  R . As  measures  of  scintillation  depth  we  use  the 
normalized  standard  deviation  of  these  quantities  defined  as 

s2  - [e(r2)  - (e(r))2]1/2/e(r) 


and 


s4  - [e(r4)  - (e(r2))2]1/2/e(r2) 

These  scintillation  indexes  are  both  used  throughout  this  report;  is 

generally  used  in  the  theoretical  developments  while  S2  is  generally 
used  In  discussions  of  measurements.  For  practical  purposes,  the  two 
indexes  are  related  to  each  other  as  S2  “ 0.5S4  . The  exact  relation- 
ship depends  upon  the  statistical  distribution  of  R . This  relationship 
for  some  representative  distributions  and  some  alternative  scintillation 
indexes  are  discussed  In  the  appendix. 
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For  the  case  of  very  large  z0(20  "*■  *) 
lowing  results: 


Mercier 
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derives  the  fol- 


E(R  ) - A A* 
o o 


(tru°  *or  all  z ) 
o 


E(R4)  - AoA*[2  - «p(-2*J)| 
s4  ' 1 - 


P 2 (u,v) 
R 


exp[2^^p^(ttVv)|  ~ 1 

K)  - 1 


exp 


(3.4) 


(3.5) 

(3.6) 


(3.7) 


Mercier  shows  that,  as  zq  becomes  large,  all  of  the  even-order  moments 
of  R coincide  with  the  even-order  moments  of  the  Rice  distribution, 
which  is  the  distribution  of  the  amplitude  of  a constant  signal  plus  in- 
phase  and  quadrature  Gaussian  noise.  Mercier  then  asserts  that  R Indeed 
has  a Rice  distribution.  The  many  small  signal  components  scattered  from 
irregularities  in  the  phase  screen  appear  to  add  together  randomly  and 
Independently  (and  thus  becoming  normal  by  the  Central  Limit  Theorem) , 
while  the  unscattered  component  of  the  field  represents  the  constant  term 
to  which  the  normal  random  components  are  added  to  yield  a Rice  distribu- 
tion. Some  of  the  other  possible  distributions  for  R are  discussed  in 
the  appendix. 
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In  the  general  case,  the  expression  for  E(R  ) , and  consequently 
, as  a function  of  <t>Q  and  zq  is  a complicated  expression  which  is 
amenable  to  numerical  evaluation  only  for  _<  10  . However,  in  certain 
circumstances  approximate  expressions  can  be  developed. 


'-'T.A'T^'R-gV' 


For  small  $o  Mercler  developed  the  approximation 

s4  ■ *♦&  + -MmA’)-1  <s-s> 

which  is  valid  for  all  z . As  z approaches  infinity  we  have 

o o 

S4  « (3.9) 

2 

while  for  z «r/X  we  have 
o o 

S4  ' l7i  (Uo/roK/’  (3‘10) 

For  (j>o  > 1 , however,  a peak  occurs  in  in  the  vicinity  of 

XzQ/r^  * tt / *^3<Jio  as  a result  of  a focusing  effect  which  is  not  predicted 
by  the  approximate  expressions.  This  focusing  can  be  visualized  as  the 
result  of  the  irregularities  refracting  rays  which  converge  together  into 
one  zone.  In  the  focusing  zone,  the  scintillation  index  can  exceed  1. 

This  result  seems  odd  in  comparison  to  a Rayleigh  distribution,  in  which 
all  of  the  power  is  in  the  fluctuation,  but  for  which  - 1 . However, 
in  the  focusing  zone  the  in-phase  and  quadrature  terms  are  not  independent 
so  that  the  fluctuations  are  more  likely  to  add  in  phase  with  the  specular 
component  and  give  a larger  amplitude.  As  zq  becomes  larger,  the  com- 
ponents become  more  random  and  the  scintillation  Index  approaches  1. 

I 

Figure  3.4  illustrates  the  behavior  of  with  zq  for  selected 

values  of  <J>q  . The  focusing  effect  yielding  > 1 is  evident  for 

<f>  > 1 . This  curve  was  computed  from  the  approximate  expressions  given 

° 12 

by  Singleton,  who  used  some  results  from  geometric  optics.  Figure  3.5, 
taken  from  Singleton,  illustrates  contours  of  constant  scintillation 
index  S^  . Region  I is  the  "near  zone"  where  S^  - ; Region  II  is 

the  "far  zone"  where  S^  - fl  - exp (“2$^  ; Region  III  is  a "saturation 

zone"  where  S^  - 1 ; and  Region  IV  is  the  "focusing  zone"  where 

Xz  / r^  - ir//3<j>  and  where  S,  > 1. 
o o o 4 
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3.16 
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The  calculations  have  assumed  an  Isotropic  screen;  that  Is,  $ ■ 1 
in  Eq.  3.2  so  that  contours  of  equal  correlation  in  the  phase  screen  are 
circles.  In  general,  Briggs  and  Parkin's  theory  admits  8 > 1 , or  ellip- 
tical contours  of  equal  correlation.  For  this  case  Briggs  and  Parkin 
replace  the  approximation  of  Eq.  3.8  by 


/2fo|l  - (cos  ux  cos  u2)1/2  cos  j (U;L  + u2) 


1/2 


(3.11) 


2 2 2 

where  tan  u.  - 2Xz  /irr  and  tan  u.  - 2Xz  /ir8  r . The  effect  of  the 
loo  2 o o 

anisotropy  is  to  slow  the  rate  at  which  approaches  its  limiting  value 

as  the  observer  moves  away  from  the  screen. 


It  has  also  been  assumed  that  the  source  is  infinitely  far  from  the 

screen.  Corrections  for  a source  at  a distance  z from  the  screen  can 

s 

be  made  by  replacing  zq  in  the  preceding  expressions  by 


z'  m z z /(z  - z ) 

o o s v o s 


(3.12) 


This  result  was  derived  by  Bowhill. 
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3.3  PHASE  AND  ANGLE  SCINTILLATION 

In  addition  to  the  amplitude  scintillations  in  the  diffracted  wave, 
there  are  phase  and  angle-of-arrival  scintillations  as  well.  Actually, 
the  angle-of-arrival  scintillation  is  a consequence  of  the  rate  of  change 
of  signal  phase  with  observer  position.  The  deviation  in  ar.gle  of  arrival 
is  given  by 


6 “ (X/27t)  to  (3.13) 

where  (x)  is  the  phase  deviation  from  mean  at  the  point  of  observation. 
If  the  covariance  function  for  the  phase  deviations  is  given  by 
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$(x)  - E [4> (a.  + x)<j>(a)]  > where  <Js(x)  is  assumed  to  be  stationary,  then 
the  variance  of  angle  of  arrival  is  given  by 

E(62)  - (— X2/ 4tt2)  (3.14) 

dx 

provided  the  second  derivative  exists  at  x * 0 . In  the  general  two- 
dimensional  case  the  orthogonal  components  of  the  angle  variation  can  be 
obtained  in  the  obvious  way. 

Immediately  beyond  the  phase  screen,  all  of  the  power  in  the  fluctua- 
tions is  concentrated  solely  in  the  phase  variations  since  the  amplitude 
variations  develop  only  as  the  wave  propagates  beyond  the  screen.  Bowhill^ 
shows  that  for  the  shallow  phase  screen,  as  the  observation  plane  is  far 
from  the  phase  screen: 

• Correlation  between  phase  in  observation  plane  and  phase  at 
screen  -*■  0 

• Correlation  between  amplitude  and  phase  in  the  observation 
plane  -*■  0 

• For  the  stationary  Gaussian  covariance  of  phase  at  the  screen, 
E[<j>(x,y,zo)<Kx  + u,y  + v,zq)]  -*•  0.5E[  (x,y,0H(x  + u,y  + v,0)]; 
that  is,  the  mean  square  phase  variation  in  the  observa- 
tion screen  tends  to  one-half  that  in  the  phase  screen 

and  the  covariance  function  take  the  same  form. 

Mercier's^  expression  (3.7)  for  the  autocorrelation  function  of  ampli- 
tude for  the  shallow  phase  screen  shows  that  it  also  takes  the  same  form 
as  the  autocorrelation  function  of  phase  at  the  screen. 

3.4  NATURAL  SCINTILLATION 

Estimates  of  the  natural  scintillation  enable  us  to  view  in  proper 

perspective  the  scintillation  induced  by  the  modified  ionosphere.  Fremouw 
13 

and  Rino  of  Stanford  Research  Institute  have  recently  presented  a model 
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for  the  average  natural  F-layer  scintillation  at  VHF/UHF  which  supplies 
parameters  for  the  Briggs  and  Parkin  model.  The  model  assumes  constant 
values  for  the  parameters: 

a ■ axial  ratio  * 10 
Ah  ■ layer  thickness  ■ 100  km 
zq  ■ height  of  layer  ■ 350  km 

For  the  parameters  r (■  transverse  scale  size)  and  AN(-  RMS  fluctuation) 

formulas  are  given  which  depend  upon  geomagnetic  latitude  X in  degrees, 

8 

local  time  of  day  in  hours  t , day  of  the  year  D , and  sunspot  number 
R . The  scale  size  tq  Is  assumed  to  depend  on  latitude  only,  assuming 
the  value  300  meters  near  the  equator  and  900  meters  in  latitudes  from 
about  15  to  60  degrees.  In  the  middle  latitudes  AN  depends  upon  only 
latitude  and  time  according  to 

AN  - (6  x 108)(1  + 0.4  cos  nt/12)  exp^-(Xg  - 32.5)2/10oj  (3.15) 

where  AN  is  in  units  of  electrons  per  cubic  meter. 

For  the  geomagnetic  latitude  of  Boulder  (about  40  degrees),  we  have 

calculated  the  expected  average  scintillation  index  as  a function  of 

scan  angle  i off  vertical  in  a magnetic  north-south  plane  using  the 

values  from  the  Fremouw  and  Rino  model,  the  Briggs  and  Parkin  model  for 
2 

4>o  , and  expression  (3.11)  for  the  scintillation  index.  The  wavelength 
is  2 meters,  and  the  source  is  assumed  to  be  infinitely  far  away.  Since 
the  dip  angle  is  about  68  degrees  near  Boulder,  the  angle  to  the  mag- 
netic field  goes  to  zero  at  i = -22  degrees  . The  average  scintillation 
index  is  plotted  in  Fig.  3.6  for  local  times  of  midnight,  noon,  and 

6 a.ra.  and  p.m.  Recall  that  the  scintillation  S,  is  about  twice  the 

10  4 

index  S£  reported  by  Bcwhill. 


In  geomagnetic  equatorial  regions  the  Fremouw  and  Rino  model  gives 


AN  = (5.5  x 109)(1  + 0.05R)  |l  - 0.4  cos  tt(D  + 10)/91.25] 

• jexp |-(t/4)2 j + exp |-(t  - 23.5)2/3.52 
• jexp [-(Xg/12)2] | (3.16) 

This  expression  assumes  its  largest  values  at  the  equator  (X  * 0),  at 

© 

night,  around  the  autumnal  and  vernal  equinoxes  (D  - 80  and  266),  and 
increases  with  sunspot  number  (R  - 100  corresponding  roughly  to  a sunspot 
maximum) . For  t^OjD^SOjR^SO,  and  X ■ 0 , we  obtain 
AN  - 2.7  x io  . In  comparison,  the  model's  peak  value  for  Boulder  is 

Q 

about  4.8  x 10  , or  about  50  times  smaller. 


For  vertical  incidence  at  the  geomagnetic  equator,  the  Briggs  and 
Parkin  equation  for  $ using  the  values  from  the  Fremouw  and  F.ino  model 
at  150  MHz  (X  ■ 2 meters)  yields  <J>q  = 1.1  . Since  a ground  observer  is 
in  the  far  zone,  the  scintillations  are  fully  developed  and  one  would 

-4 

expect  - 1 . Under  the  same  conditions  at  noon,  however,  <j>Q  - 10 
so  that  the  scintillations  would  be  negligibly  small. 

Large  natural  scintillations  in  communication  satellite  signals  at 

4 and  6 GHz  are  observed  occasionally  in  equatorial  regions.  These  were 
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reported  by  Craft  and  Westerlund  and  by  Taur.  The  scintillations 
sometimes  exceed  6 to  8 dB  peak-to-peak,  or,  equivalently,  ^0.25  . 

Taur  reports  that  these  scintillations  occur  less  than  1 percent  of  the 
time,  that  their  occurrence  coincides  very  closely  with  sunset  in  the  F- 
region,  about  2000  local  time,  and  that  they  occur  most  frequently  around 
the  equinoxes. 

Models  of  average  VHF  scintillation  do  not  predict  these  strong 
scintillations  at  much  higher  frequencies.  For  example,  the  Fremouw  and 
Rino  model,  under  conditions  which  give  K 1.1  at  150  MHz  with  resul- 
tant strong  fading,  would  predict  4>q  - 0.03  at  6 GHz.  Because  of  the 
small  wavelength  at  6 GHz  and  the  large  structure  sizes,  the  ground  observer 

would  be  in  the  near  zone  and  resultant  scintillation  index  S,  would  be 

4 

about  0.005.  The  much  stronger  scintillations  observed  occasionally  at 
SHF  imply  that  the  electron  density  fluctuations  are  probably  much 
larger  than  predicted  by  the  model  and  that  the  structure  sizes  are  small 
enough  to  bring  a ground  observer  into  the  focusing  zone  or  far  field. 
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SATELLITE  TRANSMISSION  EXPERIMENTS 


The  satellite  transmission  experiments  in  PRAIRIE  SMOKE  I-V  were 

16~"22 

conducted  and  reported  by  Aeronomy  Corporation.  Transmissions  from 

a geostationary  satellite  and  from  a number  of  polar-orbiting  satellites 
were  recorded  at  ground  receiving  stations  located  so  that  the  line-of- 
sight  penetrated  the  heated  region  of  the  ionosphere.  The  geostationary 
satellite  was  ATS-5,  transmitting  at  its  telemetry  frequency  (136. A7  MHz). 
The  orbital  satellites  were  US  Navy  navigation  satellites  in  circular 
polar  orbits  at  about  1025  km  altitude.  These  navigation  satellites 
transmit  at  about  150  and  400  MHz.  In  PRAIRIE  SMOKE  III  some  measurements 
were  made  of  pulsed  signals  at  30  and  50  MHz  from  the  Lockheed  satellite 
STP-71-2,  a part  of  the  ERIS  experiment. 

Two  mobile  receiving  stations  were  equipped  with  VHF  converters, 

HF  receivers  with  linear  detectors,  and  magnetic  tape  and  strip  chart 
recorders.  One  receiving  station,  installed  in  a travel  trailer,  was 
used  m;  ily  in  the  geostationary  experiments  since  it  was  less  easily 
moved.  The  other  station,  installed  in  a Winnebago  vehicle  shell,  was 
easily  moved  and  was  used  in  the  orbital  experiments. 

Usually,  three  antennas  feeding  separate  receivers  were  deployed 
in  a triangular  pattern,  with  maximum  separation  up  to  300  meters.  The 
amplitude  of  the  signal  in  each  receiver  channel  was  recorded  on  a strip 
chart  or  magnetic  tape  or  both.  Data  on  magnetic  tape  was  later  digitized 
for  computer  processing. 

The  geostationary  experiments  were  conducted  from  receiving  sites 
located  near  Newcastle,  Lance  Creek,  and  Lusk  in  Eastern  Wyoming  (see 
map,  Fig.  4.1).  These  sites  were  chosen  to  view  the  ATS-5  satellite 
through  the  heated  region  above  the  transmitter  near  Platteville.  The 
elevation  angle  to  the  satellite  from  these  sites  is  about  40  degrees, 
so  that  the  lines  of  sight  made  about  a 30-degree  angle  with  the  magnetic 
field  in  the  heated  region. 


The  orbital  experiments  were  conducted  at  a number  of  sites  in  Wyoming , 
Colorado,  and  Nebraska  (see  map,  Fig.  4.1).  The  sites  were  chosen  so  that 
the  llne-of-sight  to  the  satellite  would  sweep  through  the  heated  region 
during  the  satellite  pass.  Many  of  the  measurements  were  made  near  the 
downfield  point  from  the  heated  region — near  Hillsdale,  Wyoming — so  that 
the  line-of-sight  to  the  satellite  could  make  the  smallest  possible  angle 
to  the  magnetic  field.  However,  because  of  the  earth's  rotation,  the 
roughly  north-south  ground  track  of  the  satellite  varied  from  pass  to 
pass  to  the  east  or  west  of  the  ideal  track  which  would  yield  a magnetic 
field  angle  going  to  zero  as  the  line-of-sight  passed  through  the  heated 
region.  Consequently,  the  smallest  field  angle  reported  is  about  3 degrees, 
but  each  satellite  pass  produced  a range  of  field  angles  of  Interest. 

4.1  STRENGTH  OF  SCINTILLATIONS 

Quite  a broad  range  of  fading  intensity  was  measured  in  the  geosta- 
tionary and  orbital  experiments.  Scintillations  were  generally  weaker 
for  the  geostationary  experiments  since  the  line-of-sight  was  well  away 
from  alignment  with  the  magnetic  field,  but  occasional  strong  scintilla- 
tions were  observed  in  the  evening.  Figure  4.2  shows  a scatter  plot  of 
scintillation  index  S2  observed  in  the  geostationary  experiments.  There 
is  considerable  variation  in  the  observed  values  due  to  changes  in  heater 
frequency  and  power  on  separate  occasions.  The  strong  scintillations 
occurring  in  the  evening  are  probably  due  at  least  in  part  to  the  natural 
changes  in  the  ionosphere  after  sunset.  In  one  series  of  measurements, 
the  natural  daytime  scintillation  was  estimated  to  be  about  0.01  with 
about  a 9-second  period.  After  heating  the  scintillation  index  was  in 
the  range  0.04  to  0.10  and  the  fading  had  about  a 6-second  period. 

Some  of  the  orbital  experiments  yielded  values  of  S2  in  the  range 
0.2  co  0.J  for  viewing  angles  near  the  magnetic  field  lines.  The  orbital 
measurements  will  be  discussed  in  more  detail  in  a later  section  in  which 
we  compare  the  experimental  results  to  the  theoretical  predictions. 
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Some  variation  of  scintillation  intensity  on  heater  power  and  heater 
frequency  was  observed;  however,  the  dependence  is  not  very  well  deter- 
mined. For  fu  (the  heater  frequency)  greater  than  1.02f  F-  (the  cri- 
tical  frequency),  very  little  modification  is  produced;  for  Q.95fQF2 
I < fu  < 1.02f  F_  , the  results  are  quite  variable,  with  occasionally  very 

^ ^ n O Z 

| strong  scintillation  produced;  and  for  0.75foF2  <.  f <_  0.95fQF2  , the 

\ heating  effect  is  approximately  constant.  The  dependence  on  power  was 

| very  difficult  to  assess,  but  the  scintillation  tends  to  increase  with 

I increased  heater  power. 

| 4.2  IRREGULARITY  SIZE,  ORIENTATION,  AND  VELOCITY 

l With  records  of  signals  received  at  the  displaced  antennas  it  was 

possible  to  estimate  the  size,  shape,  and  velocity  of  the  irregularities 
in  the  heated  region.  Since  the  signal  source  on  the  geostationary  satel- 
lite is  fixed,  the  fluctuations  in  the  recorded  signals  were  principally 
due  to  the  drift  of  irregularities  in  the  ionosphere  relative  to  the  line- 
of-sight.  By  measuring  the  average  time  between  signal  peaks  at  each 
site,  one  can  estimate  the  ratio  of  the  structure  size  to  its  velocity. 
(For  example,  a 10-second  fade  period  could  be  caused  by  a 10-meter  struc- 
ture driftinj  at  1 meter  per  second  or  by  a 1-kilometer  structure  drift- 
ing at  100  meters  per  second.)  Since  the  signals  received  at  the  dis- 
placed antennas  were  almost  perfectly  correlated  except  for  a time  shift, 
it  was  possible  to  estimate  the  drift  velocity  as  the  ratio  of  the  antenna 
separation  to  the  time  shift  between  records.  The  structure  size  could 
then  be  deduced  by  multiplying  the  size-to-velocity  ratio  by  the  estimated 
velocity.  Finally,  by  assuming  a two-dimensional  Gaussian  correlation 
function  for  the  phase  deviation  in  the  ionosphere,  one  can  then  relate 
the  deduced  structure  size  to  the  dimensions  of  the  correlation  ellipse. 

In  contrast  to  the  geostationary  measurements,  the  fluctuations  in 
the  received  signal  for  the  orbital  satellites  were  due  mainly  to  the 
motion  oi  the  line-of-sight  relative  to  the  irregularities  as  the  satei  • 
lite  moved  along  its  orbit  at  about  7200  m/s.  The  natural  structure  drift 


rate  of  about  50  m/s  was  negligible  in  comparison.  Since  the  satellite 
orbit  was  known  very  precisely,  the  rate  of  motion  of  the  line-of-sight 
through  the  heated  region  could  be  computed  and  the  else  and  orientation 
of  the  irregularities  deduced  from  the  measured  period  of  the  fluctua- 
tions. Since  the  line-of-sight  changed  from  pass  to  pass,  and  since  the 
measurement  sites  were  moved  from  near  the  downfield  point  to  several 
other  sites,  line-of-sight  sections  of  the  heated  region  from  several 
different  aspect  angles  were  observed.  Consequently,  the  size  and  field 
alignment  of  the  irregularities  and  the  size  and  shape  of  the  whole  region 
could  be  estimated.  A histogram  of  the  structure  sizes  observed  in 
PRAIRIE  SMOKE  II  and  IV  is  shown  in  Fig.  4.3. 

The  conclusions  regarding  irregularity  size  and  velocity  are  sum- 
marized below: 

• About  80  percent  of  the  time,  the  transverse  structure  size 
rQ  lay  in  the  range  75  to  300  meters,  averaging  about  100 
meters. 

• Some  fine  structure  of  the  order  of  10  meters  was  observed, 
confirming  some  results  obtained  in  radar  backscatter 
measurements . 

• Drift  velocity  is  of  the  order  of  20  to  50  m/s  in  an  east- 
west  direction. 

• The  irregularities  are  elongated  along  the  magnetic  field 
line,  having  a slight  banana  shape  tilting  slightly  more 
vertical  at  the  top  and  more  horizontal  at  the  bottom. 

4.3  COMPARISON  OF  EXPERIMENTS  WITH  THEORY 

In  order  to  compare  the  experimental  results  to  the  theoretical 

predictions,  it  is  convenient  to  recast  the  Briggs  and  Parkin  formula 

(Eq.  3.1)  into  a slightly  different  form.  The  plasma  frequency  f is 

2 2 w 

related  to  the  electron  density  by  f^  * Nc  r e h , where  N is  the  elec- 
tron density,  c is  the  speed  of  light,  and  rg  is  the  classical  electron 
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Figure  4.3.  Structure  Sizes  Observed  in  PRAIRIE  SMOKE  II  and  IV 
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radius.  If  we  equate  the  critical  frequency  t0  fc^e  plasma  frequency 

at  maximum  electron  density  N and  recall  that  c ■ fX  , we  obtain  by 

' m 

substitution  into  3.1  the  result 


if,  - 

l' 


r. 

E 

S" 

n 


„ 5/2 /f  fA4  E(AN2)Ahar  sec  i 

,2  TT  I O Z|  O 
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(4.1) 


This  form  is  more,  convenient  in  that  the  absolute  value  of  RMS  electron 
density  deviation  can  be  regarded  more  simply  as  a fraction  of  the  maxi- 
mum electron  density. 

In  order  to  get  a feeling  for  the  numbers  in  Eq.  4.1,  let  us  substi- 
tute the  following  values: 

A * 2m 
f - 150  MHz 

f F,  = 7.5  MHz 
o z 

Ah  = 100  km 
a = 10 

r = 100  m 
o 


Then 


2 2.73  x 103  E |(AN/Nin)^|sec  i 

° (100  sin2  i|>  + cos2  ty)*^ 


(4.2) 


* 4 
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For  viewing  along  the  field  line,  $ « 0 and  i **  22  degrees  , so  that 
2 3 2 

d>  - 2.94  x 10  E[(AN/N  ) ] ; for  vertical  viewing,  4>  *»  22  degrees  and 

° 2 m 2 2 
i « 0 so  that  q>  - 7.07  x 10  E[(AN/N  ) ] . The  fractional  deviation 
o m 

in  electron  density  to  achieve  scintillation  index  S,  - 0.25  is  about 

-3  -3 

3.3  x 10  for  viewing  along  the  field  and  6.6  x 10  for  vertical 


About  the  only  means  to  compare  the  predictions  of  Briggs  and  Parkin 
theory  with  the  experimental  measurements  is  in  the  variation  of  scintil- 
lation index  with  the  angle  ip  between  the  line-of-sight  and  the  magnetic 
field.  Since  only  ip  and  the  zenith  angle  vary  with  line-of-sight  in 
the  Briggs  and  Parkin  formula,  and  since  in  the  far  field  the  scintilla- 


tion index  is  directly  proportional  to  $ , we  have 
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S,  « sec*' “ ilsin  ’J)  + — =■  cos  ip 


-1/4 


(4.3) 
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Also,  since  the  variation  of  sac*'  * i within  about  10  degrees  of  the 
field  line  is  less  than  10  percent,  we  can  neglect  its  effect  and  have 
instead 


« yjin^'  ip  + -^2  cos^  ipj  (4.4) 

Clearly  this  function  has  a peak  at  ip  *=  0 for  the  axial  ratio 

a > 1 . The  rate  of  falloff  from  this  peak  depends  upon  a , being  very 

1/2 

slow  for  small  a and  approaching  falloff  as  esc  ip  for  large  values 
of  a . Figure  4.4  illustrates  this  variation. 

For  several  of  the  orbital  satellite  experiments  the  scintillation 

index  S2  was  plotted  as  a function  of  both  time  and  field  angle  ip  . 

Results  from  one  such  experiment  are  replotted  in  Fig.  4.5.  Because  of 

the  offset  of  the  satellite  tracks  from  the  heated  region  the  field  angle 

ip  was  never  smaller  than  about  3 degrees . The  heavy  dashed  lines  across 

1/2 

Fig.  4.5  represents  falloff  with  ^ as  esc  ip  . Although  that  is  the 
fastest  rate  of  falloff  from  the  peak — corresponding  to  an  infinite  axial 
ratio — the  measured  falloff  is  somewhat  taster.  The  Briggs  and  Parkin 
prediction  does  not  fit  the  measured  variation  with  angle  very  well. 
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5 A SCINTILLATION  MODEL  FOR  THE  MODIFIED  IONOSPHERE 

One  of  the  assumptions  of  the  Briggs  and  Parkin  scintillation  model 
is  that  the  layer  of  the  ionosphere  producing  the  scintillation  is  hori- 
zontally infinite  with  stationary  statistics.  However,  the  heated  region 
of  the  ionosphere  has  a horizontal  extent  comparable  to  the  diameter  of 
the  heater  beam — about  80  km  at  half  power  at  300  km  altitude.  Also, 
since  the  degree  of  modification  depends  on  the  intensity  of  the  heater 
beam,  we  should  expect  maximum  effect  at  the  center  of  the  heated  region 
with  a decline  in  both  horizontal  and  vertical  directions  away  from  the 
center. 

A consequence  of  the  finite  size  and  shape  of  the  heated  region  is 
that  the  phase  deviation  along  a line-of-sight  should  depend  not  only  on 
the  zenith  angle  1 and  the  field  angle  \p  according  to  the  Briggs  and 
Parkin  formula,  but  also  on  how  close  the  line-of-sight  passes  to  the  cen- 
ter of  the  heated  region.  To  account  for  this  additional  variation,  we 
have  modified  the  Briggs  and  Parkin  scintillation  model  to  account  for 
both  horizontal  and  vertical  variation  of  electron  density  fluctuations 
in  the  heated  region.  (This  development  follows  the  suggestions  outlined 
by  Drs.  S.A.  Bowhill  and  E.K.  Walton  of  Aeronomy  Corporation  in  a private 
communication.) 

5.1  A MODIFIED  SCINTILLATION  FORMULA 

2 

Suppose  that  the  variance  E(AN  ) of  the  electron  density  fluctua- 
tions in  the  slab  of  ionosphere  is  given  by 

E(AN2)  - f (x,y,z)  ; zmin  < z < zmax  (5.1) 

where  the  characteristic  dimensions  of  its  variation  are  assumed  to  be 

large  relative  to  the  dimensions  r and  ar  of  the  striations.  The 

o o 

geometry  is  indicated  in  Fig.  5.1,  with  an  observer  at  the  origin  of  the 
coordinate  system  viewing  a source  along  a line-with-direction  cosines 
(xg/rg,  yg/rs,  zg/rg)  where  rg  - xg  + yg  + zg  . The  magnetic  field  B 
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Figure  5.1.  Coordinate  System  for  Modified  Scintillation  Model 


has  a dip  angle  6 and  lies  in  the  yz-plane,  so  that  the  y-axis  points 
toward  magnetic  north.  The  line-of-sight  to  the  source  makes  an  angle 
to  the  magnetic  field  and  a zenith  angle  i ■ cos  ^(Zg/rg)  to  the 
vertical. 

At  each  point  along  the  line-of-sight,  parameterized  by  altitude  h 
as  (x,y,z)  - (hXg/Zg,  hys/Zg,  h)  we  assume  that  there  exists  a horizontal 
slab  of  incremental  thickness  dh  (corresponding  to  the  Ah  in  Briggs 
and  Parkin* s expression)  contributing  an  incremental  random  phase  shift 
d<p  . We  assume  that  the  mean  square  variation  in  the  electron  density  in 
the  slab  is  given  by  the  value  E(AN2)  - f(x,y,z)  - f(hXg/Zg,  hyg/zg,  h) 
appropriate  for  the  point  at  which  the  line-of-sight  pierces  the  slab. 

Then  to  obtain  the  total  mean  square  phase  variation  contributed  along  the 
line-of-sight  we  integrate  to  obtain 

z 

max 

$2  = E(4>2)  = 7r1//2r2A2aroB-1  sec  i J'  f(hXg/zg,  hyg/Zg,  h)  dh  (5.2) 

zmin 

This  reduces  exactly  to  the  Briggs  and  Parkin  formula  (3.1)  for  the  case 

f = constant  and  Ah  - z - z , . 

mdx  min 

We  shall  assume  that  the  correlation  function  in  the  horizontal  plane 
for  the  emerging  wave  is  the  stationary  expression 

P^(?»n)  - exp(-52/r2  - n2/ <0rQ) 2) 

as  before,  ignoring  the  horizontal  variation  in  f(x,y,z). 


Let  ua  assume  for  the  function  f the  form 


This  peculiar-appearing  form  attempts  to  match  the  falloff  in  electron 
density  variance  due  to  the  effects  of  the  ionospheric  heater  beamshape. 
Because  the  heater  beam  is  circular,  f has  circular  symmetry  in  each 
horizontal  plane  z » constant  . However,  because  the  "blobs"  of  electron 
density  become  extended  along  the  magnetic  field,  the  locus  of  centers 
of  symmetry  for  various  values  of  z is  a line  parallel  to  the  magnetic 
field.  At  each  altitude  z , the  center  of  symmetry  is  displaced  (xm»ym) 
from  the  point  of  intersection  of  the  magnetic  field  line  passing  through 
the  origin  with  the  constant  altitude  plane.  The  scale  distance  for  1/e 
falloff  in  each  horizontal  plane  is  R ; the  scale  distance  for  vertical 
1/e  falloff  is  H . The  contour  of  constant  variance — say  l/e — is  an 
ellipsoid  with  circular  horizontal  sections.  Figure  5.2  illustrates  the 

1 . ellipsoid  for  H <<  R in  the  observer's  coordinate  system.  In  gene- 
ral the  ellipsoid  could  be  chopped  off  at  the  ends  by  the  z * zmax  or 

z = zm^n  planes.  For  the  case  H ->  ~ the  ellipsoid  degenerates  to  a 

cylinder  having  a circular  horizontal  cross  section,  skewed  along  the 

field  lines  and  truncated  by  the  z » z . and  z = z „ planes. 

min  max 

Given  the  form  (5.3)  for  f(x,y,z)  , the  integral  in  Eq.  (5.2)  can 
Vo  vivalv  • . The  result  is 
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where 


A ■ 

[xs  + <ys  * zsctn6)2|/(zsR2) 

+ 1/H2 

(5.5) 

B - 

■2j[XmXS  + ym(yS  + zSctn6)|/ 

(2sr2)  + 2„/h2| 

(5.6) 

C - 

/ 2 U.  2\/D2  2 /u2 

(x  + y )/R  + z /H 
\ m m 

(5.7) 

This  expression  can  be  shown  to  reduce  to  the  Briggs  and  Parkin  formula 
in  the  case  R -*■  » and  H -*•  ® , or  equivalently,  as  the  horizontal  and 
vertical  variation  is  eliminated. 

? 2 

Equation  5.4  can  be  recast  using  the  relation  (f  F~>"  ~ N c r h 

o L m e 

as  before  to  yield 


3 

*o  * ^2  <f0Vf)4“ro8‘1(sec 


x A exp(-C  + B2/4A) 


X jcr£lAl/2<W  + B/2A>]  - erflAl/2(2„in  + B/2A)I| 


(5.8) 


where  ®(A^ax)  eclua^  to  ^ in  Eq.  5.4  or  to  the  peak  mean  square 
deviation  at  the  center  of  the  heated  region. 


It  is  frequently  convenient  to  estimate  the  maximum  phase  deviation 
that  could  be  introduced  along  a line-of-sight  through  the  disturbed 
region.  That  maximum  occurs  for  a line-of-sight  parallel  to  the  magnetic 
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field  line  directly  through  the  center  of  the  region.  In  that  case  we 


xm  “ ym 
xg  - 0 


-zg  ctn  3 


4>  * 0 


0 ■ a 
A - 1/H2 


B - -2z  /R 
m 

C - z2/H2 
m 


lx  we  assume 


2 ->  -oo 

min 


we  have: 


ffoF2r 


r H<sec  i)E 
o 


(5.9) 


5.2  COMPARISON  WITH  BRIGGS  AND  PARKIN  FORMULA 

Some  example  cases  were  computed  to  compare  the  scintillation  pre- 
dicted by  the  modified  theory  with  that  predicted  by  the  Briggs  and  Parkin 
formula.  Figure  5.3  illustrates  the  angular  dependence  in  the  magnetic 
N-S  plane  for  three  cases:  R * °°  , H ■ « (no  variation,  Briggs  and 
Parkin  case);  R ■ 50  km  , H * <*>  (no  vertical  variation);  R ■ H ■ 50  km 
(both  horizontal  and  vertical  variation).  The  center  of  the  heated  region 
was  assumed  to  lie  directly  up  the  magnetic  field  line  from  the  observer 
(that  is,  xm  * * 0) . The  falloff  with  angle  for  the  finite  region  is 

clearly  much  more  rapid  than  for  the  Briggs  and  Parkin  case.  Also,  the 
peak  scintillation  index  is  smaller  when  vertical  variation  in  the  elec- 
tron density  fluctuation  is  included.  Figure  5.4  illustrates  the  varied 
rate  of  falloff  from  the  peak  for  various  values  of  the  1/e  radius  R ; 


SCINTILLATION  INDEX  (LINEAR  SCALE) 
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the  vertical  variation  is  suppressed  (H  ■ »)  since  it  mainly  affects  the 
peak  value  and  not  the  shape  of  the  curve. 

The  effect  of  moving  the  center  of  the  heated  region  away  from  the 

direct  upfield  point  is  illustrated  in  Fig.  5.5.  The  peak  scintillation 

index  is  depressed  and  the  peak  occurs  at  a higher  angle  for  northward 

shift  y . A shift  of  y kilometers  means  that  the  observer  would  have 
m m 

to  move  northward  y^  kilometers  to  view  the  center  of  the  heated  region 
directly  along  a magnetic  field  line. 

5.3  PARAMETERS  FOR  THE  MODEL 

Now  that  we  have  a model  (Eq.  5.8)  for  computing  the  phase  variance 
along  any  line-of-sight  as  well  as  rules  for  relating  that  variance  to  the 
scintillation  index  (Eq.  3.11),  we  lack  only  means  for  determining  the 
appropriate  parameter  values  to  use  in  the  model. 

The  Aeronomy  Corporation  experiments  have  provided  estimated  values 
for  the  scale  size  rQ  and  axial  ratio  a . Figure  4.3  illustrates  the 
observed  range  of  values  for  rQ  , with  the  most  frequently  observed  value 
being  about  100  meters.  This  is  considerably  smaller  than  the  range  from 
dOO  to  900  meters  for  the  transverse  scale  size  in  the  natural  ionosphere. 
The  axial  ratio  a is  about  10  for  either  the  artificially  modified  iono- 
sphere or  the  natural  ionosphere. 

The  signal  frequency  f and  the  corresponding  wavelength  X are 

presumed  to  be  known.  The  critical  frequency  varies  continually,  usually 

in  the  range  5 to  10  MHz.  For  simple  calculations,  we  can  assume  that 

f F0  « 7.5  MHz. 
o i 


The  remaining  parameters  of  the  modal  describe  the  variance  in  elec- 
tron density  in  the  region,  the  size  and  shape  of  the  region,  and  its 
location  relative  to  the  observer. 
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A simple  model  for  predicting  the  dimensions  and  peak  electron  den- 

23 

sity  deviation  in  the  heated  region  was  given  by  Pendyala  and  Thome  on 

the  basis  of  their  RF  backscatter  measurements.  According  to  their  model, 

the  peak  electron  density  variance  occurs  at  the  center  of  the  heater  beam 

at  an  interaction  altitude  z for  which  the  heater  frequency  equals  the 

tn 

local  plasma  frequency.  (Typically  this  interaction  altitude  is  of  the 

order  of  225  to  275  km.)  The  peak  relative  electron  density  variance 

e(an2  /N2)  is  given  by 
\ max  m/ 

e(aN*  /N*)  = 106P2/9  (5.10) 

\ max  m/ 

2 

where  P is  the  peak  power  density  (in  W/m  ) in  the  heater  beam  at  the 

interaction  altitude.  On  the  basis  of  their  RF  scattering  measurements, 

-4 

Pendyala  and  Thome  reported  a saturation  level  of  about  10  corresponding 
-5  2 

to  P *>  3 x 10  W/m  . No  similar  saturation  effect  was  observed  in  the 
scintillation  measurements.  However,  agreement  with  the  square  law  power 
dependence  of  E(AN^ax/Nm)  ’ or  e<luivalent:ly » linear  variation  of  scintil- 
lation index  with  power,  was  reported  by  Aeronomy  Corporation  in  PRAIRIE 
20 

SMOKE  IV.  For  purposes  of  our  model,  we  shall  ignore  the  saturation 

-4 

effect,  but  values  much  larger  than  10  will  be  viewed  somewhat 
suspiciously. 

The  vertical  extent  of  the  heated  region  is  assumed  to  be  very 
small,  with  about  15  kilometers  between  vertical  1/e  altitudes;  this 
corresponds  to  H * 7.5  km  in  Eq.  5.3.  With  such  rapid  vertical  varia- 
tion, the  exact  values  for  z . and  z are  not  very  important.  We 
’ min  max 

can  assume  values  of  200  and  350  km,  respectively;  these  values  effectively 
correspond  to  -«  and  +»  , and  these  latter  values  could  be  used  in  the 
formulas  with  negligible  difference. 

If  we  assume  a Gaussian  shape  for  the  heater  main  beam,  we  can  com- 
pute the  horizontal  scale  distance  R for  the  heated  region.  Suppose 
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the  relative  power  pattern  for  the  heater  beam  is  given  by 
P(0)  a exp (-0^/0^)  , where  0 is  the  angle  off  vertical  and  0g  is  the 
angle  at  which  the  power  is  down  to  1/e  of  the  peak  value  (0g  - 0.6  x 3 dB 
beamwidth) . In  the  plane  of  the  interaction  altitude  z , the  angle  off 


vertical  is  8 


\ ^ (Vx2  + y2/'z  j + y2/z  for  small  angles. 


Then  in  the  interaction  altitude  plane  we  have 


P « exp  - 


2 . 2, 


Since  the  electron  density  variance  is  proportional  to  P , we  find 


E(AN2  /N2)  « exp {- 
\ max  m / M 


« exp f-(x2  + y2)/R2 


so  that 


R = z 0 //2 
m e 


(5.11) 


The  two  remaining  parameters  x^  and  ym  define  the  horizontal 

shift  of  the  center  of  the  heated  region  from  the  geomagnetic  field  line 

passing  through  the  observation  point.  The  heated  region  lies  above  the 

heater  transmitter,  but  the  center  may  be  shifted  north  or  south  slightly 

by  the  bending  of  the  0-  or  X-raode  transmissions.  However,  since  we  are 

generally  interested  in  cases  where  the  observer  is  directly  downfield 

from  the  center  of  the  heated  region,  we  can  assume  x = y =0  in  most 

m m 

cases  of  interest. 

5. A COMPARISON  OF  THE  MODIFIED  THEORY  WITH  EXPERIMENTS 

Let  us  now  compare  the  computed  results  from  the  refined  scintilla- 
tion model  with  some  of  the  measured  scintillation  data  from  the  satellite 
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experiments  in  PRAIRIE  SMOKE  III.  The  measurements  were  taken  on  three 
satellite  passes  on  8 September  1972.  The  amplitude  of  the  satellite's 
150-MHz  signal  was  recorded  on  each  pass  at  the  Pine  Bluffs,  Wyoming, 
receiver  site  operated  by  Aeronomy  Corporation.  Concurrent  measurements 
on  the  last  two  passes  were  also  made  at  the  Hillsdale  site  a few  miles 
away.  The  amplitude  scintillation  index  S£  was  then  computed  at  time 
steps  of  about  5 seconds.  These  values,  reported  in  Ref.  18,  are  replctted 
here  in  Figs.  5.6,  5.7,  and  5.8  for  the  three  satellite  passes. 

The  angle  i|/  between  the  line-of-sight  to  the  satellite  and  the 
geomagnetic  field  ranged  from  a minimum  of  about  3 to  10  degrees  up  to 
angles  exceeding  90  degrees.  Only  the  values  less  than  about  15  degrees 
correspond  to  significant  scintillation  levels. 

The  parameter  values  used  for  the  calculations  are  summarized  in 
Table  5.1.  Some  of  the  assumed  values  were  chosen  as  discussed  in  the 
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Figure  5.6.  Comparison  of  Model  with  Experimental  Data,  Satellite  Pass  1 
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Figure  5.7,  Comparison  of  Model  with  Experimental  Data,  Satellite  Pass  2 
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Figure  5.8.  Comparison  of  Model  with  Experimental  Data,  Satellite  Pass  3 


TABLE  5.1 


PARAMETER  VALUES  FOR  MODEL  CALCULATIONS 
FOR  PRAIRIE  SMOKE  III  EXPERIMENTS 


Parameter 

Value 

Parameter 

Value 

(Assumed  or 
Calculated) 

ro 

100  m 

0 

a 

10 

Zm 

275  km 

R 

35  km 

EKax/Nm) 

4.8  x io"4 

H 

7.5  km 

(Known  or  Measured) 

Zmin 

200  km 

f 

150  MHz 

Z 

max 

350  km 

X 

2 m 

X 

m 

0 

foF2 

5.5  MHz — Pass 
7.8  MHz — Pass 

7.4  MHz— Pass  3 


preceding  section.  R and  E (AN^gv/N^ ) were  computed  from  the  known 
parameters  of  the  heater  transmitter — 1.9  MW  power,  15  dB  gain,  and  a 
16-degree  half-power  beamwidth — using  Eqs.  5.10  and  5.11.  The  remaining 
parameters  were  either  known  beforehand  or  measured  for  each  pass. 

The  calculations  were  made  for  lines-of-sight  in  the  magnetic  north- 
south  vertical  plane  passing  through  the  center  of  the  heated  region. 

The  observer  was  assumed  to  lie  directly  downfield  from  the  center  of  the 
region.  For  this  case  the  field  angle  \p  goes  to  zero,  in  contrast  to 
the  actual  measurements  where  the  minimum  was  a few  degrees.  The  calcula- 
tions yield  slightly  different  values  for  Sj  along  lines-of-sight  north 
and  south  of  the  field  line  because  of  the  sec  i variation  in  the  phase 
variance  (Eq.  5.8). 


In  addition  to  the  idealized  geometry,  calculations  were  made  which 
attempted  to  duplicate  the  actual  geometry  for  the  satellite  pass.  The 
calculations  yielded  minimum  angles  within  about  1 degree  of  those  reported 
for  the  measurements.  (Actually,  the  field  angles  reported  in  Ref.  18 
were  computed,  too,  and  not  measured.)  However,  the  scintillation  indexes 
were  negligibly  different  from  the  idealized  geometry. 


The  agreement  between  the  computed  values  and  the  measurements  is 
good  in  Figs.  5.6  and  5.7,  but  rather  poor  in  Fig.  5.8.  The  agreement  is 
best  for  the  cases  where  the  scintillation  is  strong  and  the  field  angle 
\p  is  small.  These  cases  are  clearly  the  ones  of  greatest  interest  in 
any  possible  systems  applications.  Because  of  the  few  cases  available  for 
validating  the  model,  it  should  be  used  with  some  caution. 
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SCINTILLATION  EFFECTS  ON  RADAR  SYSTEMS 


We  have  seen  in  the  preceding  sections  that  artificial  heating  of 
the  ionosphere  can  produce  modest  levels  of  scintillation  along  lines-of- 
sight  passing  through  the  heated  region.  In  addition  to  its  use  as  an 
experimental  mechanism  in  studying  the  ionosphere,  artificial  modifica- 
tion of  the  ionosphere  could  potentially  be  used  to  degrade  the  perfor- 
mance of  RF  systems  requiring  one-way  or  two-way  lines  of  sight  through 
the  F-region.  The  two  such  systems  that  come  immediately  to  mind  are 
long-range  radars  for  tracking  missiles  or  satellites  in  space  and  satel- 
lite-to-ground  (or  vice  versa)  communication  relays. 

The  objectives  of  this  section,  which  treats  the  radar,  and  the  next 
section,  which  treats  satellite  communication,  are  to  identify  the  degra- 
dations produced  by  ionospheric  scintillation  and  to  relate  that  degrada- 
tion parametrically  to  the  scintillation  index.  By  treating  the  scintil- 
lation level  parametrically,  we  can  extend  many  of  the  results  to  natural 
or  nuclear-induced  scintillation. 

Two  principal  effects  on  radar  operation  are  produced  by  scintilla- 
tion. The  first  is  a change  in  the  detection  statistics,  generally  a 
degradation  of  high  detection  probabilities  and  an  enhancement  of  low 
detection  probabilities.  In  Sec.  6.2  the  detection  probabilities  for 
various  signal-to-noise  ratios  are  computed  as  a function  of  the  scintil- 
lation index  assuming  a Rice  distribution  for  the  scintillation  amplitude. 

The  implications  of  that  assumption  are  discussed  in  Sec.  6.1.  The  second 
effect  is  a degradation  of  angular  tracking  accuracy  due  to  fluctuations 
in  the  angle  of  arrival  of  the  radar  signal.  Tracking  errors  are  evaluated 
using  a Monte  Carlo  computer  simulation  for  tracking.  Some  results  from 
this  simulation  are  presented  in  Sec.  6.3. 

6.1  THE  AMPLITUDE  DISTRIBUTION  FUNCTION 

There  is  some  dispute  concerning  the  correct  distribution  function 
to  use  in  describing  the  amplitude  scintillations  which  occur  in  the  pre- 
sence of  striations.  Some  of  the  alternatives  are  discussed  in  the  appendix. 
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Mcrcler  has  asserted,  and  other  researchers  (e.g.,  Briggs  and  Parkin, 

12 

Singleton  ) have  agreed,  that  the  distribution  In  the  limit  of  large  dis- 
tances from  the  phase  screen  Is  Rician.  Mercler's  result,  although  the 
derivation  may  not  be  mathematically  rigorous,  Is  intuitively  satisfying 
as  was  pointed-out  in  Sec.  3.2.  On  the  other  hand,  some  later  workers 
(see  for  example  Rino  and  Fremouw  ) have  been  dissatisfied  with  this 
assumption.  They  have  found  that  the  measured  data  for  the  normal  iono- 
sphere are  not  well  fitted  by  a Rician  distribution,  but  rather  require 
a more  complicated  function  based  on  two-dimensional  Gaussian  statistics 
in  which  the  in-phase  and  quadrature  components  of  the  scattered  power 
have  unequal  sigmas.  See  the  appendix  for  a further  discussion  of  this 
problem.  Rino  and  Fremouw  found  that  their  data  implied  a transverse 
scale  size  r^  - 600  meters  , and  inferred  that  the  far-field  conditions 
under  which  the  Rician  amplitude  distribution  of  Mercier  was  derived  would 
rarely  be  appropriate  at  typical  radar  operating  frequencies. 

For  the  purposes  of  the  present  discussion,  however,  it  is  impor- 
tant to  note  that  the  artificially  induced  strlatlons  produced  by  iono- 
spheric heating  appear  to  have  significantly  smaller  transverse  scale 
sizes  than  those  found  in  the  natural  ionosphere.  A glance  at  Fig.  4.3 
reveals  that  the  distribution  of  transverse  scale  sizes  measured  in  the 
PRAIRIE  SMOKE  II  and  IV  experiments  peaks  at  around  rQ  * 100  meters  , 
that  50  percent  of  the  time  rQ  < 150  meters  , and  that  90  percent  of  the 

time  r < 300  meters  . Now,  Fig.  3.5  shows  that  far-field  conditions 
° 2 

may  be  assumed  to  apply  whenever  XZQ/ro  > 2 • Assuming  that  the  phase 
3creen  is  at  least  300  km  from  the  radar  and  that  the  target  is  far  be- 
yond the  screen,  it  is  easy  to  compute  the  maximum  value  of  rQ  for 
which  the  far-field  approximations  are  generally  valid,  as  a function  of 
the  radar  frequency  f . This  has  been  plotted  in  Fig.  6.1,  on  the  follow- 
ing page.  It  is  obvious  from  this  figure  that  Rino  and  Fremouw  were  cor- 
rect in  asserting  that  a ground-based  receiver  looking  at  a target  near 
the  zenith  will  generally  not  be  in  the  far-field  with  respect  to  natural 
ionospheric  striations  with  a mean  transverse  scale  size  of  the  order  of 
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600  meters.  On  the  other  hand,  it  is  equally  obvious  that  a radar  operat- 
ing at  VHF  or  UHF,  and  even  at  L-band  or  up  to  S-band,  will  generally  be 
in  the  far-field  with  respect  to  the  artificially  induced  striations. 

Accordingly,  for  the  present  application  it  is  not  unreasonable  to 
accept  Mercier’s  hypothesis  and  use  his  Rician  amplitude  probability- 
density  function.  For  a one-way  propagation  path  through  the  phase-screen, 
the  amplitude  density  function  is 


53 


exp  -(1  + e)(R/R 


2 /g(l  + 8)  R/R 


(6.1) 


where  R is  the  observed  amplitude,  Rq  is  the  value  which  would  have 

been  observed  in  the  absence  of  scintillations,  and  6 is  a parameter 

which  is  described  below.  The  function  I (x)  is  the  modified  Bessel 

o 

function  of  zero  order. 


The  parameter  8 which  appears  in  Eq.  6.1  is  related  to  the  funda- 
mental properties  of  the  phase  screen.  In  Sec.  3.1  we  defined  the  phase 
function  <J>(x,y)  at  the  screen,  and  its  RMS  value 


^E(<|>2) 


We  now  define  the  parameter  8 to  be 


(6.2) 


1 - e 


It  can  be  shown  (R.W.  Hendricx  3)  that  the  power  in  the  undiffracted 
("specular")  portion  of  the  signal  transmitted  through  the  phase  screen 
is  proportional  to  exp(-<£2^  , and  the  powpr  in  the  diffracted  ("random") 
portion  of  the  signal  is  proportional  to  1 - exp|-<|>2j  . The  parameter 
8 therefore  is  equal  to  the  ratio  of  the  specular  to  random  power  in  the 
received  signal;  it  may  be  thought  of  as  a kind  of  signal-to-noise  ratio. 
Alternatively,  it  is  a measure  of  the  strength  of  the  phase  screen:  a 

large  value  of  8 corresponds  to  weak  diffraction  and  a small  value  of 
8 corresponds  to  strong  diffraction. 


It  is  not  surprising,  therefore,  that  8 is  related  to  the  scintil- 
lation indices  defined  in  Sec.  3.2  and  the  appendix.  From  Eq.  3.6  it  is 
easily  seen  that 

s«  ■ 1 - (-nb)2  <6-3> 

Thereupon  the  known  relationship  between  and  S£  for  the  Rice  dis- 
tribution (see  Fig.  A.l)  can  be  used  to  compute  6 as  a function  of  S£  . 
Figure  6.2  shows  and  3 as  functions  of  the  scintillation  index 


For  radar  problems  in  which  the  signal  is  propagated  twice  through 

the  phase  screen,  a slightly  different  amplitude  distribution  function, 

P2(R;Rq,8)  , must  be  used.  This  is  found  from  p^  by  making  use  of  the 

fact  that  the  two-way  distribution  of  R is  equivalent  to  the  one-way 
2 

distribution  of  R . Standard  techniques  then  permit  derivation  of  the 
function  P2  : 


p2(R;Ro’6)  = 6Xp 


-(1  + 0)  - 0 


j 2^0(1  + B)R/R 


(6.4) 


Here  again,  R is  the  observed  instantaneous  value  of  the  ampli- 
tude, and  Rq  is  the  value  which  would  have  been  observed  in  the  absence 
of  striations.  The  parameter  6 , which  characterizes  the  phase  screen 
rather  than  the  process,  retains  the  same  value  whether  for  one-way  or 
two-way  propagation.  It  should  be  noted  here  that  the  amplitude  R is 
related  to  the  measured  signal-to-noise  ratio  (S/N)  in  either  case  by 
the  identity 


R2  = 2 (S/N) 


(6.5) 
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6 and  RMS  Phase  Angle  Fluctuation 
tion  Index  S£ 


A H-426S4 


6.2  RADAR  DETECTION  PROBABILITY  AND  DETECTION  RANGE 

26 

Standard  radar  detection  theory  (see  for  example  DiFranco  and  Rubin, 
especially  Chapter  9)  computes  the  detection  probability  for  typical  signal 
processing  in  the  following  manner.  Assume  a steady  sinusoidal  signal  of 
known  amplitude  and  unknown  phase,  added  to  zero-mean  Gaussian  noise  of 
known  power.  The  radar  receiver  performs  an  envelope  detection  and  then 
determines  the  presence  or  absence  of  the  signal  by  comparing  the  peak 
voltage  of  the  video  signal  with  a pre-selected  threshold.  The  probabi- 
lity that  the  threshold  will  be  exceeded  when  the  signal  is  indeed  pre- 
sent defines  the  detection  probability,  P^. 

The  threshold  level,  r^  , is  selected  to  provide  a desired  proba- 
bility of  false  alarm,  P^  ; this  is  simply  the  probability  that  the  ran- 
dom noise  alone,  in  the  absence  of  the  signal,  will  produce  a voltage 
peak  which  exceeds  the  threshold.  It  is  of  course  desirable  to  minimize 
the  false  alarm  probability,  but  the  price  that  must  be  paid  for  a very 
low  value  of  is  an  increased  possibility  that  'iven  when  the  signal 

is  present  a threshold  crossing  will  fail  to  occurs  in  other  words,  a 
reduced  probability  of  detection.  In  practical  radar  systems  a compro- 
mise must  be  made  whereby  an  acceptable  value  of  P-  is  chosen;  values 

-4  -8  1 

between  10  and  10  are  commonly  used. 

It  can  be  shown  that  the  probability  density  that  the  peak  video 
amplitude  will  be  R when  the  amplitude  of  the  sinusoidal  signal  is  RQ  , 
is  given  by  the  Rician  distribution 

p(R,R0)  - R exp | - (r2  + R^)/2jlo(RoR)  (6.6) 

Here  R and  R are  normalized  to  the  noise  voltage,  so  R is  given 
o 

by  £q.  6.5  and 

R2  = 2 (S/N)  (6.7) 

o o 
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where  (S/N) o is  the  known  signal-to-noise  ratio  for  the  steady  input 
sinusoidal  signal.  This  represents  the  case  of  a target  whose  cross  sec- 
tion is  intrinsically  non-fluctuating. 


For  a given  false  alarm  probability, 
must  be  defined  by  requiring  that 


the  amplitude  threshold  level 


t 

f- 

\ 

1 


| which  leads  at  once  to  the  expression 

\ r,  = f-2  InP, 

t b t 

P 


(6.8) 


E 


Having  set  the 
becomes 


threshold  accordingly,  the  detection  probability  then 


f 


P 


d 


f p(R,Rq)  dR 


(6.9) 


If  the  sinusoidal  signal  amplitude  is  not  constant,  but  fluctuates 
either  because  the  target  cross-section  fluctuates  or  because  the  propa- 
gation medium  induces  scintillations  (or  both) , then  the  detection  proba- 
bility must  be  computed  by  multiplying  Eq.  6.9  by  the  amplitude  probabi- 
lity density  distribution  f(RQ)  and  integrating  overall  values  of  Rq 
i * e . , 


00  00 

pd •//  p<r-r;>£(r;>  dR dRi 

° rb 


(6.10 
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We  consider  the  case  of  a non-fluctuating  target  but  a propagation 
medium  which  contains  a phase  screen  consisting  of  a layer  of  striations 
that  causes  amplitude  scintillations.  For  a two-way  path  through  this 
medium  the  amplitude  scintillations  have  a distribution  function  of  the 
form  of  p2  as  given  by  Eq.  6.4.  Thus 


00  00 

Pd(R0,8)  - If  p(R,RpP2(R^Ro,6)  dR  dR^  (6.11) 

0 V 

This  expression  has  been  evaluated  as  a function  of  scintillation 
index  S2  (through  the  parameter  0)  and  vacuum  signal-to-noise  ratio 
(S/N)q  (through  the  parameter  Rq)  for  several  typical  values  of  the  false 

alarm  probability  P,  . Some  of  these  results  are  shown  in  Figs.  6.3  and 

1 -4  -8 

6.4:  Fig.  6.3  for  Pf  - 10  and  Fig.  6.4  for  Pf  = 10  . Notice  that 

the  scintillations  tend  to  increase  the  detection  probability  for  targets 

of  low  signal-to-noise  ratio  but  to  decrease  the  detection  probability 

when  the  signal-to-noise  ratio  is  large.  A similar  effect  has  long  been 

observed  for  cases  of  fluctuating  targets  in  a neutral  environment. 

These  results  can  be  used  to  compute  degradation  factors  which  indi- 
cate how  the  scintillations  reduce  the  range  at  which  a radar  can  achieve 
a given  detection  probability  on  any  target.  The  degradation  factors  can 
be  expressed  in  a form  which  is  independent  of  the  target  or  its  signal- 
to-noise  ratio,  but  which  depends  only  on  S2»  , and  P^. 

We  start  with  the  standard  radar  equation  in  the  form 

S/N  ■*  (S/N)1(R1/R)4a/a1  (6.12) 

where  here 


S/N  is  the  measured  signal-to-noise  ratio 
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(S/N)^  is  the  minimum  signal-to-noise  ratio  required  for  de- 
tection by  the  given  radar 

is  the  range  at  which  a target  of  cross-section 
will  just  produce  this  minimum  signal-to-noise  ratio 

* 

R is  the  current  range  to  the  target, 
a is  the  current  cross-section  of  the  target 

The  parameters  (S/N)^,  R1  , and  o^  are  given  or  may  be  computed  from 
other  intrinsic  parameters  defining  a specific  radar. 

New,  in  the  absence  of  striations  ("vacuum")  it  is  possible  to  com- 
pute the  minimum  signal-to-noise  ratio,  (S/N)q  , which  is  required  to 
assure  a specified  detection  probability,  P^  , with  a given  false  alarm 
probability  . Curves  like  those  in  Figs.  6.3  and  6.4  can  be  used  to 

find  (S/N)  . Then,  for  a target  of  known  cross  section,  the  radar  equa- 

o £ 

tion  can  be  used  to  compute  the  maximum  detection  range,  Rq(0)  * at 
which  this  target  can  be  expected  to  be  detected  with  probability  P^  . 

It  is  given  by 


R (0) 
o 


»>1  q 
(S/N)o  Oj 


When  striations  are  present  (S^  > 0) , the  signal-to-noise  ratio  re- 
quired to  achieve  the  same  detection  probability  will  be  changed  to 
(S/N)^  which,  as  we  have  seen,  may  be  larger  or  smaller  than  (S/N)q  . 
Similarly,  the  detection  range  will  be  changed  to 


VS2> 


*1; 
(S/B) ' o 


1 


£ 

There  should  be  no  danger  of  confusing  this  R and  Rq  with  the  symbols 
R and  R which  represent  signal  amplitudes  in  other  expressions. 
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Clearly,  therefore,  the  ratio  of  the  detection  ranges  with  and  without 
the  scintillations  is 


w 

Ro(°) 


(S/N) 

(s7n) 


11/4 


o 

* 

0 


(6.13) 


and  is  independent  of  the  specific  parameters  which  describe  the  radar 
or  the  target. 

In  this  manner  we  have  computed  the  degradation  factor,  the  ratio 

R (S_)/R  (0)  , as  a function  of  S0  for  selected  values  of  P,  and  for 
° 2 o _/  ^ o u 

three  values  of  : 10  ,10  , and  10  . The  results  are  shown 

plotted  in  Figs.  6.5,  6.6,  and  6.7,  respectively. 

It  is  clear  from  these  results  that  a radar  can  suffer  severe  degra- 
dation of  its  single-pulse  detection  range  on  a non-fluctuating  target 
when  the  intervening  medium  has  a scintillation  index  S2  greater  than 
about  0.3.  The  results  remain  the  same  for  multiple-pulse  detection  pro- 
vided that  the  integration  period  is  less  than  the  scintillation  correla- 
tion time.  For  fluctuating  targets,  or  for  multi-pulse  detection  when 
the  integration  period  is  greater  than  the  scintillation  correlation  time, 
one  would  expect  in  general  that  the  degradation  would  be  less  significant 
than  that  found  above,  but  the  specific  results  would  depend  on  the  sta- 
tistics for  the  case  of  a fluctuating  target. 


6 . 3 RADAR  TRACKING 

When  a radar  tracks  a target  which  passes  behind  a striated  region, 
the  scintillations  cause  tracking  errors  which  increase  the  prediction 
error  at  the  intercept  point.  Although  the  details  are  somewhat  depen- 
dent on  the  specific  scenario,  it  is  possible  to  express  the  results  in 
a form  which  is  relatively  scenario  insensitive,  and  to  parameterize  some 
of  the  intrinsic  variables.  This  is  what  we  shall  do  in  the  following 
pages. 
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Figure  6.5.  Ratio  of  Detection  Range  to  Vacuum  Detection  Range  as  Function  of 

-4 

Scintillation  Index,  for  Non-Fluctuating  Target  and  ■ 10 
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The  scintillations  which  affect  the  tracking  results  are  of  two 
kinds:  amplitude  scintillations  such  as  we  have  already  discussed,  and 

angle  scintillations  which  will  be  described  below.  These  cannot  be  sepa- 
rated in  the  analysis  because  the  amplitude  scintillations  cause  changes 
in  the  natural  radar  errors  in  range  and  angle  in  addition  to  their  effect 
on  detection  range. 


6.3.1  Radar  Errors  as  Influenced  by  Amplitude  Scintillations 

A radar  usually  measures  three  coordinates  of  the  target  at  each 
pulse:  the  range  (R) , and  two  angular  coordinates,  the  azimuth  (A)  and 

the  elevation  (E) , or  sine-space  coordinates  (u,v).  In  any  measurement 
of  R,A,E  , the  radar  will  make  unavoidable  errors.  In  addition  to  bias 
errors,  which  by  definition  are  constant  and  unknown  during  any  given 
track,  there  are  stochastic  errors  whose  statistics  at  least  are  known. 

The  statistical  parameters  can  usually  be  expressed  in  a form  which  exhi- 
bits the  fact  that  they  have  a part  which  is  independent  of  the  signal- 
to-noise  ratio  and  another  part  which  is  signal-to-noise  dependent.  Thus, 
we  write  that  any  measurement  of  the  parameter  X (where  X may  be  R, 

A,  or  E)  has  a zero-mean  Gaussian  error  of  standard  deviation  a , where 


2 

a 

x 


Ts7nT 


(6.14) 


The  six  coefficients  (^x*By)  are  known  from  the  calibration  measurements 
performed  on  any  given  radar.  Approximate  but  reasonably  accurate  values 
for  can  also  be  computed  from  the  known  parameters  of  the  radar 

antenna  and  waveform. 


When  tracking  a target  in  "vacuum,'1  the  signal-to-noise  ratio  is 
given  by  the  radar  equation  as  a function  of  the  target  and  radar  param- 
eters and  the  range  to  the  target  at  any  given  time.  Thus  the  variances 
2 

ax  can  be  estimated  and  used  in  the  tracking  logic,  but  the  random  mea- 
surement errors  for  each  pulse  are  not  known;  this  leads  to  inescapable 
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tracking  errors  and  an  eventual  target  prediction  error  whose  statistics 
are  often  estimated  by  performing  a Monte-Carlo  simulation  of  the  track- 
ing process. 

When  strlations  are  present  along  the  radar  line-of-sight,  the  value 
of  (S/N)  is  itself  a random  variable  whose  statistics  are  given  by  Eqs. 

6.4  through  6.6.  In  this  case  the  only  way  to  ascertain  the  eventual  tar- 
get prediction  error  is  to  perform  the  Monte-Carlo  calculation  using 

\ 

samples  drawn  randomly  from  the  appropriate  distributions.  The  algorithm 
by  which  this  has  been  done  in  the  present  investigation  will  be  discussed 
in  Sec.  6.3.3  below. 

6.3.2  Angular  Scintillations 

In  addition  to  the  "normal"  angular  errors  described  by  Eq.  6.14, 
the  presence  of  an  Intervening  striated  region  also  Induces  fluctuations 
in  the  angle-of-arrival  of  the  radar  returns  as  measured  at  the  receiver. 
Representing  the  striated  region  as  a phase  screen  again,  and  assuming 
that  far-field  conditions  are  applicable,  the  angular  scintillations  may 
be  described  by  a zero-mean  Gaussian  distr  jution  about  the  "vacuum" 
angle,  with  variance 

a]  = 2(X/2Trro)2^  (6.15) 

By  making  use  of  Eq.  3.6  and  the  fact  that  S^  = 2S£  (see  the  appen- 
dix) , we  can  write 

$2  « 2S ^ [ 1 + 2S 2 1 

an  approximation  which  upon  comparison  with  Fig.  6.2  is  seen  to  be  very 
good  on  the  range  0 <_  S£  ^ 0.45  and  adequate  out  to  S£  - 0.5  . Sub- 
stituting this  into  Eq.  6.15  leads  to  the  expression 
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zyi  + 2s; 

V(ro/X) 


(6.16) 


Since  the  striations  are  assumed  to  be  elongated  along  the  magnetic 
field  lines,  with  a "longitudinal*'  scale  size  much  greater  than  rQ  , it 
must  be  recognized  that  the  angular  scintillations  characterized  by 
occur  only  in  the  direction  normal  to  both  the  radar  line-of-sight  and 
the  magnetic  field  lines  in  the  region  of  the  phase  screen.  Any  scintil- 
lations parallel  to  the  field  lines  would  be  of  negligible  magnitude, 
and  will  be  ignored  in  the  present  analysis.  Accordingly,  the  angular 
scintillations  are  characterized  by  their  standard  deviation,  , and 
their  direction,  u , in  the  radar  viewing  plane  normal  to  the  line-of- 
sight.  (Numerical  values  of  a are  given  later  in  Fig.  6.11.) 


6.3.3  The  Monte-Carlo  Tracking  Simulation 

In  order  to  evaluate  the  effects  of  both  amplitude  and  angle  scin- 
tillations on  the  radar  tracking  function,  a Monte-Carlo  tracking  simula- 
tion program  was  written.  This  program  accepted  a given  target  trajectory 
and  specified  radar  parameters.  It  placed  a heated  patch  (striated  region) 
of  appropriate  size,  shape,  and  altitude  at  a specified  location  with 
respect  to  the  trajectory,  so  as  to  interfere  with  part  of  the  tracking; 
it  was  geometrically  impossible  to  interfere  with  the  entire  track  using 
only  one  patch. 


The  program  simulated  the  tracking  process  by  using  a random  number 
generator  in  conjunction  with  the  known  statistical  properties  where  re- 
quired to  represent  the  noise,  radar  errors,  and  scintillation  effects. 

In  order  to  avoid  extraneous  issues,  a perfect  tracking  filter  was  assumed. 
After  making  tracking  measurements  for  a specified  time  period  at  a given 
tracking  rate,  using  the  data  corrupted  by  the  effects  previously  described, 
the  program  projected  the  calculated  trajectory  forward  in  time  to  a spe- 
cified time  or  altitude  (as  desired)  and  estimated  the  target  location. 


This  was  compared  with  the  true  target  position  at  this  time  to  yield  the 
prediction  error,  D . The  above  tracking  and  projection  procedure  was 
repeated  a given  number  of  times,  using  different  random  numbers  at  each 
iteration,  and  the  results  were  accumulated  to  produce  the  mean  error, 

D , and  its  standard  deviation,  . The  entire  procedure  could  be  car- 
ried out  for  a range  of  input  values  of  t'.  i parameters  S2,  r^/X  . and  p. 

A flow  chart  for  this  program  is  shown  in  Fig.  6.8. 

6.3.4  The  Baseline  Parameters 

In  using  the  Monte-Carlo  Tracking  Simulation  program,  specific  sce- 
narios had  to  be  used.  There  was  no  way  to  avoid  choosing  particular  tra- 
jectories and  radar  parameters  for  each  calculation.  However,  the  results 
have  been  expressed  in  a relative  form  which  is  intended  to  reduce  the 
specific  scenario  dependence;  in  addition  several  of  the  most  significant 
parameters  were  separately  varied  over  their  natural  range  to  evaluate 
their  individual  effects  on  the  results.  It  is  hoped  that  this  form  of 
analysis  will  provide  a picture  of  the  tracking  effects  due  to  the  scintil 
lations  which  is  relatively  independent  of  scenario.  This  subject  will 
receive  additional  comment  in  the  discussion  to  follow. 

It  was  Intended  to  run  a series  of  different  trajectories  and  radar 
types  through  the  program.  Unfortunately,  the  program  turned  out  to  have 
such  a long  running  time  that  it  was  impractical  to  carry  out  this  proce- 
dure. We  were  able  to  use  only  one  of  the  planned  trajectories  and  radar 
locations,  although  we  were  able  to  study  the  effects  of  changing  the 
radar  parameters  by  varying  some  of  them  in  separate  calculations. 

The  trajectory  covered  an  8,000-km  ground  range,  with  the  radar 
location  as  its  intended  impact  point.  Figure  6.9  illustrates  the  nature 
of  this  trajectory  relative  to  the  radar  by  showing  the  target  range,  ele- 
vation angle,  and  altitude  (above  the  ground)  as  functions  of  the  time 
after  launch  in  seconds . 
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Figure  6.8,  Flow  Chart  for  Monte-Carlo  Tracking  Simulations 


mpmm  miu  mm^ti*u^ZYV 


/S£3M (VV 


xj 

s 


o 

4J 

(U 

*-» 

3 

(S 

>* 

M 

O 

V 

o 

V 

*cd 

t-i 

H 

& 


«wi  ‘(a)  nm 


f 


The  radar  used  In  the  baseline  calculations  was  assumed  to  operate 
at  150  MHz,  the  frequency  at  which  the  measured  scintillation  data  were 
obtained.  As  a practical  matter,  radars  used  to  track  targets  in  space 
typically  operate  at  somewhat  higher  frequencies,  often  in  the  UHF  band, 
but  the  parameters  of  the  assumed  radar  are  representative  of  this  type 
of  radar  except  for  frequency.  These  parameters  are  given  in  Table  6.1. 

This  radar  has  a detection  range  (in  vacuum)  on  a l-square  meter 
target  of  = 7,900  km  . Its  effective  measurement  error  sigmas  (see 
Sec.  6.3.1)  are: 

= 12  m , B„  = 150  m 

K 

Angle:  Afl  = 0.6  mrad  , B = 15.0  mrad 


Range : A^ 


TABLE  6.1 
RADAR  PARAMETERS 


Frequency 
Wavelength 
Peak  Power 

Average  Power 

Pulse  Length 

Pulse  Compression  Ratio 

Effective  Range  Resolution 
Antenna  Beam  Width 
Antenna  Gain 
Detection  Threshold 

System  Noise  Temperature 

System  Losses 
Horizon  Limit 


f 

150  MHz 

X 

2 m 

p 

o 

5 MW 

P 

250  kW 

T 

200  us 

F 

pc 

200 

6R 

300  m 

69 

1.4  deg 

G 

44  dB 

(S/N)o 

10  dB 

tn 

600K 

L 

6 dB 

E 

0 

2 deg 
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A heated  patch  (striated  region)  with  typical  dimensions  was  placed 
at  the  appropriate  altitude  in  two  different  positions  relative  to  the 
radar  and  target  trajectory.  These  two  cases,  shown  in  Fig.  6.10,  illus- 
trate the  two  extreme  conditions  for  artificially  induced  scintillations 
placed  to  interfere  with  radar  tracking  of  a specific  trajectory  by  a spe- 
cific radar;  they  are: 

• Case  1,  Scintillations  First.  The  patch  is  placed  so  as  to 
obscure  the  beginning  of  track  (at  E = Eq) , but  the  last  por- 
tion of  track  will  be  in  the  clear. 

• Case  2,  Scintillations  Last.  The  patch  is  placed  so  as  to 
obscure  the  end  of  track,  but  the  first  portion  will  be  in 
the  clear. 

In  each  case,  the  times  (in  seconds  relative  to  launch  time)  at  which 
scintillations  would  begin  and  end  for  the  given  trajectory  and  radar 
location  were  precomputed  and  used  as  entries  in  the  simulation  program. 

In  all  calculations  the  tracking  rate  was  held  constant  at  1 pulse 
per  second,  and  the  correlation  time  for  the  striations  was  quasi- 
arbitrarily  set  at  5 seconds. 

Four  different  values  of  the  relative  scale  size  r /X  were  used: 

o 

rQ/X  = 200,  100,  50,  25  . Since  this,  is  the  only  place  where  the  radar 
wavelength  actually  enters  the  calculations,  this  parameterization  has 
the  effect  of  removing  any  significance  from  the  fact  that  the  baseline 
radar  happened  to  be  given  an  operating  frequency  of  150  MHz.  This  range 
of  values  of  rQ/X  spans  the  region  of  practical  interest.  For  rQ/X  > 200 
(high  frequency  or  large  striations)  the  effect  on  radar  tracking  is,  as 
we  shall  see,  quite  negligible.  And  rQ/X  < 25  (low  frequency  or  small 
striations)  will  practically  never  be  observed,  given  the  distribution  of 
striation  scale  sizes  shown  in  Fig.  4.3.  It  should  be  noted  that  this 
entire  range  of  relative  scale  sizes  satisfies  the  far-field  condition 
for  a heated  region  at  300  km  altitude. 


74 


To  get  an  idea  of  the  magnitudes  of  the  angle  scintillation  errors, 

we  have  computed  a from  Eq.  6.18  as  a function  of  S„  for  the  four 

o 2 

selected  values  of  relative  scale  size.  These  are  shown  plotted  in  Fig. 
6.11.  For  comparison  we  have  plotted  on  the  same  figure  some  values  for 
the  standard  deviation  of  the  natural  ("vacuum")  radar  measurement  angle 
error,  a A from  Eq.  6.14,  for  several  values  of  S/N  , using  the  param- 
eters of  the  baseline  radar.  Among  other  things,  this  figure  illustrates 

why  the  effects  of  angular  scintillations  with  r /X  > 200  are  so 

o 

negligible. 

In  the  baseline  calculation,  tracking  was  carried  out  at  one  pulse 
per  second  from  T * 960  seconds  (E  = 2 degrees)  to  T = 1245  seconds 
(E  = 11.8  degrees).  This  last  track  time  was  selected  so  that  the  tra- 
jectory could  be  projected  ahead  300  seconds  (untracked)  to  a point  at 
about  400  km  altitude.  The  results,  when  stated  in  the  relative  form  to 
be  described  below,  should  be  insensitive  to  the  specific  value  of  the 
last  track  time  (projection  time  was  one  of  the  parameters  which  were 
later  .varied) . 

A pre-computation  showed  that  the  mean  miss  distance  varied  with 
the  angle-error  direction  y in  a uniform  manner  over  a factor  of  2 as 
y went  from  0 to  90  degrees,  with  y = 45  degrees  representing  the  mid- 
way position.  Accordingly,  the  baseline  calculations  were  performed 
with  y = 45  degrees  . This  parameter  was  later  varied. 

All  calculations  were  performed  for  values  of  the  scintillation 
index  on  the  range  0 <_  S£  ^ 0.50. 

The  minimum  number  of  iterations  necessary  to  provide  consistent 
results  was  determined  by  pre-computation  to  be  80.  All  track  simulations 
were  iterated  80  times,  and  the  mean  prediction  error  D and  its  standard 
deviation  were  calculated  for  each  run. 


STANDARD  DEVIATION  OF  ANGLE  ERROR  (o.  OR  cfl),  mrad 


All  results  were  normalized  to  the  value  of  the  mean  prediction 
error,  Dq  , for  vacuum  conditions  (no  scintillations).  The  relative 
mean  prediction  error  D/Dq  is  used  to  show  the  basic  qualitative  trends 
and  approximate  quantitative  effects  due  to  variation  of  the  significant 
parameters.  In  this  relative  form  the  results  should  be  insensitive  to 
those  elements  of  the  specific  scenario  which  do  not  directly  affect  the 
variation  of  5 as  a function  of  . Insofar  as  we  were  able  to  deter- 
mine, only  the  location  of  the  heated  patch  in  relation  to  the  beginning 
and  end  of  track  has  a significant  effect  on  this  function,  and  this 
"parameter"  was  spanned  by  using  Cases  1 and  2 separately  in  all  calculations. 

6.3.5  Results  and  Parameter  Variations 

The  results  of  the  baseline  calculations  for  Cases  1 and  2 are  shown 
in  Figs.  6.12  and  6.13,  as  functions  of  scintillation  index  and  relative 
scale  size.  Notice  that  there  is  practically  no  effect  for  r^/X  = 200  . 


SCINTILLATION  INDEX, 


Figure  6.12.  Relative  Mean  Prediction  Error  as  Function  of  Scintillation  Index 
for  Several  Values  of  Scale  Size  Over  Wavelength,  Case  1: 
Scintillations  First 
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Figure  6.13.  Relative  Mean  Prediction  Error  as  Function  of  Scintillation  Index 
for  Several  Values  of  Scale  Size  Over  Wavelength,  Case  2:  . 
Scintillations  Last 


In  general,  D/Dq  increased  with  S2  , the  effect  being  negligible  for 
rQA  ■ 200  and  growing  more  significant  as  r /A  decreased.  Case  2 
(scintillations  last)  showed  impressively  greater  increase  in  the  rela- 
tive mean  prediction  error  with  S2  than  was  found  for  Case  1 (scintil- 
lation first).  This  should  be  expected,  since  in  Case  1 the  tracking 
filter  has  an  opportunity  to  recover  from  the  errors  induced  by  the  scin- 
tillations before  making  its  projection,  whereas  in  Case  2 it  has  no  such 
opportunity. 

As  a matter  of  interest,  the  value  of  5 for  the  baseline  cases 

o 

was  approximately  2 km. 

It  should  be  pointed  out  that  all  tracking  error  simulations,  how- 
ever they  are  performed,  show  quite  a lot  of  scatter.  Even  in  vacuum, 
the  standard  deviation  of  the  prediction  error  is  usually  a significant 
fraction  of  the  mean  value.  As  an  illustration  of  this  effect,  Fig.  6.14 
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SCINTILLATIONS  FIRST 

SCALF  SIZE/WAVELENGTH  (rQ/>.)  = 50 


Figure  6.14. 


SCINTILLATION  INDEX,  $2 

Relative  Mean  Prediction  Error  as  Function  of 
Scintillation  Index 


shows  the  one-sigma  error  bars  on  the  relative  prediction  error  for  one 
particular  pass  through  the  tracking  simulation.  Note,  however,  that 
the  basic  trend  is  still  clearly  apparent.  The  reader  should  imagine 
error  bars  like  these  on  all  of  the  curves  of  D/Dq  plotted  in  this 
section;  they  were  omitted  only  because  they  would  have  cluttered  the 
figures  into  incomprehensibility. 


As  a separate  calculation,  the  prediction  time  after  track  was 
varied  over  the  range  1 s <_  T <_  400  s for  one  particular  set  of  the 
other  basic  parameters  to  ascertain  the  effect  of  this  parameter  on  the 
baseline  results.  Although  the  calculation  was  only  performed  for  this 
one  case,  the  qualitative  trend  should  apply  to  all  of  the  baseline 
results,  and  there  is  no  reason  to  suppose  that  the  relative  quantita- 
tive values  should  be  much  different  for  other  cases.  Figure  6.15  shows 
the  mean  prediction  error,  relative  to  the  baseline  value  for  T = 300 
seconds  as  a function  of  projection  time.  I*-  is  quite  obvious  that  the 
mean  error  varies  linearly  with  the  projection  time. 


Figure  6.15.  Relative  Mean  Prediction  Error  as  Function  of  Prediction  Time 
After  Track 


In  another  separate  calculation,  the  angle-error  direction  u was 
varied  over  the  range  0°  <_  y £ 90°  for  one  particular  set  of  the  other 
basic  parameters  to  examine  the  effect  of  this  parameter  on  the  baseline 
results.  In  "real  life,"  the  parameter  y iB  determined  by  the  direc- 
tion of  a vector  transverse  to  the  magnetic-field  lines  in  the  heated 
patch,  projected  onto  the  plane  normal  to  the  radar  line-of-sight.  Figure 
6.16  shows  the  mean  prediction  error,  relative  to  the  baseline  value  for 
y ■ 45  degrees  , as  a function  of  angle  y . The  effect  of  this  geometri- 
cal factor  is  readily  apparent. 

Finally,  we  must  take  note  of  the  fact  that  the  radar  parameters 

might  be  quite  different  than  those  given  in  Table  6.1,  or  the  target 

2 

cross  section  might  be  different  from  the  value  a ■ 0.4  m used  in  the 
baseline  calculation,  either  because  the  target  itself  is  different  or 
because  the  same  target  is  being  observed  at  a different  radar  frequency. 
Many  of  the  radar  parameters  appear  lumped  into  the  constant  which 
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CASE  2;  SCINTILLATIONS  LAST 
S2  * 0.25 


ANGLE-ERROR  DIRECTION  WITH  RESPECT  TO  AZIMUTH/ELEVATION  AXIS,  (j. ) , deg  ' 

Figure  6.16.  Relative  Mean  Prediction  Error  as  Function  of  Angle  Error 
Direction 


is  used  in  the  condensed  form  of  the  radar  equation  given  in  Eq.  6.12. 
2 

Using  = 1 m , we  can  rewrite  that  equation  in  the  form 


S/N  » KR" 


where 


K i (S/tDjRj  §- 


Now,  by  varying  K over  an  appropriate  range  of  values  we  can  account 
for  variations  in  the  target  cross  section  and/or  many  of  the  radar 
parameters. 


Accordingly,  in  another  separate  calculation  K was  varied  over 
15  4 18  4 

the  range  10  km  £ K <_  10  km  for  one  particular  set  of  the  other 
basic  parameters.  (The  baseline  value  of  K is  K ■ 1.56  x 10^  km^.) 


The  effect  on  the  baseline  results  of  varying  this  lumped  parameter  is 

shown  in  Fig.  6.17,  where  we  have  plotted  the  relative  mean  prediction 

error,  normalized  to  its  value  for  K ■ K , as  a function  of  K . Note 

0 17 

that  the  curve  becomes  effectively  flat  for  K >_  5 * 10  , which  is  a 

good  reason  why  there  is  no  point  in  extending  the  calculation  beyond 

18 

K =*  10  . The  reason  for  this  behavior  is  easy  to  find:  for  the  power- 

ful radars  or  strong  targets  represented  by  these  large  values  of  K , 
the  radar  measurement  error  sigmas  (which  are  the  primary  beneficiaries 
of  the  increased  signal-to-noise  ratio)  are  approaching  their  limiting 
values  for  S/N  >>  1 . In  the  other  direction,  the  dashed  portion  of  the 
curve  for  K < 2 x 10^  indicates  the  region  where  the  returns  are  so 
weak  that  for  at  least  part  of  the  tracking  range  they  failed  to  cross 
the  receiver  threshold.  If  the  figure  were  extended  much  below  K “ 10^"* 
it  would  be  found  that  few  or  no  successful  track  pulses  would  be  received 
during  the  time  interval  over  which  we  are  trying  to  track,  so  that  the 
prediction  error  rapidly  approaches  infinity. 

It  is  obvious  that  the  radar  tracking  error  is  more  serious  when 
the  heated  region  obscures  the  terminal  portions  of  the  track  than  when 
it  occurs  toward  the  beginning.  In  principle,  significant  increases  in 
prediction  error  can  occur  for  such  "scintillations  last"  cases  even  when 
the  scintillation  index  is  only  moderately  severe  (0.2  < S2  0.3).  In 
practice,  however,  geographical  considerations  would  generally  make  it 
difficult  or  impossible  to  place  the  striated  patch  close  enough  to  the 
radar  to  ensure  that  no  good  track  returns  can  be  obtained  after  the 
scintillation  interval,  and  without  such  placement  the  effect  approaches 
a less  critical  "scintillations  first"  situation. 

Since  the  angle  y is  not  usually  subject  to  control,  the  results 
illustrated  by  Fig.  6.16  are  of  primarily  theoretical  interest. 

The  radar  power  vis-a-vis  target  cross  section,  on  the  other  hand, 
can  clearly  be  of  great  significance  in  evaluating  the  effecf.3  of 
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Relative  Mean  Prediction  Error  as  a Function  of  Radar  Parameter 


scintillations  on  track  performance.  We  have  seen  that  it  is  of  diminish- 
ing usefulness  to  increase  the  power  beyond  a certain  point,  but  that 
tracking  can  be  seriously  degraded  by  the  heated  patch  if  the  radar  is 
marginal  in  power  to  begin  with. 


7 SATELLITE  COMMUNICATION  SYSTEMS 

The  basic  elements  of  a satellite  communication  system  consist  of 
a ground  terminal  transmitter,  a satellite-borne  repeater,  and  a ground 
terminal  receiver.  One  or  both  of  the  "ground"  terminals  may  in  fact  be 
in  an  aircraft,  but  the  principle  is  the  same.  The  satellite  equipment 
usually  provides  only  the  functions  of  frequency  translation  and  amplifi- 
cation; signal  processing  and  demodulation  generally  occur  only  at  the 
ground  terminals. 

The  uplink  and  downlink  each  provide  a one-way  path  which  might 
intersect  a heated  patch  and  suffer  some  degradation.  It  is  geometrically 
unlikely,  though  not  impossible,  that  both  links  should  be  affected  by  a 
single  patch.  In  this  section  we  will  compute  the  degradation  for  a single 
one-way  link;  extension  to  account  for  striation  effects  on  both  links  is 
obvious . 

Each  link  is  usually  required  to  handle  multiple  simultaneous  com- 
munication channels.  Some  types  of  channels  handle  voice  messages,  others 
deal  with  digital  information.  The  common  method  used  by  satellite  sys- 
tems for  multi-channel  voice  communications  is  frequency  division  multiplex 
frequency  modulation  (FDM-FM) . To  transmit  digital  information  over 
satellite  communication  links  there  are  basically  two  commonly  used  modu- 
lation methods:  frequency-shift-keying  (FSK)  and  phase-shift-keying  (PSK) . 

Each  of  these  in  turn  may  be  classified  according  to  the  method  of  detec- 
tion used.  For  FSK  the  detection  may  be  coherent  (CFSK)  or  non-coherent 
(NFSK) . For  PSK,  one  may  again  find  coherent  detection  (CPSK)  or  non- 
coherent detection  used,  though  in  the  latter  case  what  is  usually  employed 
is  the  co-called  differentially  coherent  system  (DPSK) . For  details  con- 
cerning the  theory  and  operation  of  these  (as  well  as  other)  systems,  a 

textbook  on  communication  theory  should  be  consulted  (see  for  example 

27 

Schwartz,  Bennett,  and  Stein  ). 


The  source  of  system  degradation  which  we  shall  consider  here  is 
the  "signal  fading"  caused  by  relative  motion  of  the'  satellite  and  the 
striations  which  lie  across  the  propagation  path.  Even  in  the  absence  of 
artificially  induced  striations,  satellite  communication  links  are  sub- 
ject to  natural  fading,  and  they  are  designed  with  the  capacity  to  per- 
form adequately  under  given  expected  levels  of  fading.  Depending  on  geo- 
graphical location  and  local  time-of-day,  natural  fading  of  1 to  10  seconds 
duration  and  a depth  of  4 to  6 dB  occurs  quite  frequently,  with  longer 

and/or  much  deeper  fading  occurring  correspondingly  less  often.  Paulson 
28 

and  Hopkins  report  fading  at  250  MHz  of  the  order  of  10  dB  or  greater 
peak-to-peak  occurring  5 to  10  percent  of  the  time — mainl>  at  night — over 
a 3-year  period  at  receiving  sites  located  near  the  equator. 

Some  discussion  of  the  fading  phenomena  associated  with  the  artifi- 
cially induced  striations  may  be  found  in  the  appendix.  There  it  is  pointed 

29 

out  that  with  the  Rice  or  Nakagami  distributions  for  signal  amplitude 
almost  any  fading  depth  is  possible  with  some  small  probability.  However, 
the  fading  depth  at  any  particular  level  of  probability  can  be  related  to 
the  scintillation  index  by  comparing  the  amplitudes  at  two  specified  per- 
centile levels  of  the  distribution  function.  This  was  done  for  the  Nakagami 
distribution  and  index  in  Fig.  A. 2 of  the  appendix.  Since  it  can 

be  seen  from  Fig.  A.l  how  close  the  Nakagami  distribution  is  to  the  Rice 
distribution  on  the  one  hand  and  to  the  approximation  = 25^  on  the 

other.  Fig.  A. 2 can  be  Interpreted  in  terms  of  index  S2  , if  desired,  by 
dividing  the  given  index  in  half.  The  meaning  of  these  curves  is  that 
for  a given  value  of  the  scintillation  index  the  indicated  fading  depth 
(in  dB)  will  be  observed  5 percent  of  the  time  for  the  bottom  curve,  1 
percent  of  the  time  for  the  middle  cutve,  or  0.1  percent  of  the  time  for 
the  upper  curve. 

A perhaps  more  direct  measure  of  system  degradation  may  be  found 
by  considering  the  effect  of  the  fading  on  a digital  communication  system. 

A digital  information  signal  is  composed  of  a sequence  of  symbols,  and 
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each  symbol  Is  transmitted  as  a sequence  of  elements.  The  Probability  of 

Element  Error,  Pe  , is  defined  as  the  probability  that  the  receiver 

assumes  any  element  other  than  the  one  transmitted.  In  Chapter  7 of 

27 

Schwartz,  Bennett  & Stein,  it  is  shown  that  P is  a function  only  of 

e 

the  signal-to-noise  ratio  on  the  link,  and  expressions  for  these  functions 
are  derived  in  the  cases  of  the  four  types  of  digital  communication  sys- 
tems discussed  above.  For  DPSK  and  NFSK  systems  the  expressions  are  par- 
ticularly simple: 


For  DPSK:  Pg  = | exp[-(S/N) ] (7.1) 

For  NFSK:  Pg  - j exp  [~  ~ (S/N)]  (7.2) 


Suppose  that  the  signal-to-noise  ratio  at  a given  instant  would  have 
been  S/Nq  in  the  absence  of  striations,  but  that  because  of  fading  it 
has  the  received  value  S/N  . Using  symbols  previously  introduced,  we 
can  write 

(S/N)q  = R2/2  (6.7) 

and 


(S/N)  = R2/2  (6.5) 

The  fading  due  to  the  striations  has  the  amplitude  distribution 
Pl(R;Ro,8) 
given  in  Eq.  6.1. 
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Accordingly,  the  probability  of  element  error  in  the  presence  of 
fading  can  be  computed  by  integrating: 


PeCR0,8> 


/ 


Pe(R)Pl(R;Ro,0)  dR 


(7.3) 


For  DPSK  and  NFSK  systems  this  integral  can  be  done  analytically  and  the 
results  expressed  in  simple  closed  form.  We  find  that 


For  DPSK:  % - | i £ ^ y 


exp 


By, 


1 + B + Y, 


(7.4) 


For  NFSK:  P - $ 

e Z 


1 + 0 


1 + 8 + f yo 


exp 


0Yo/2 


1+B+f  Yo 


(7.5) 


where  y = (S/N) 

o o 

Notice  that  these  two  expressions  are  identical  except  for  a factor  of  2 
in  signal-to-noise  ratio. 

Since  the  DPSK  method  is  extremely  common  ' use,  let  us  now  focus 

onto  Eq.  7.4.  In  Fig.  7.1  we  show  P as  a ft  .on  of  scintillation 

£ 

index  S£  for  selected  values  of  the  "vacuum"  signal-to-noise  ratio  yq  , 
and  in  Fig.  7.2  we  show  Pg  as  a function  of  y for  selected  values 

of  S2  • Finally,  we  can  calculate  a link  degra'  ion  loss  factor  which 
sums  up  the  entire  situation  in  the  following  manner.  For  many  practical 
communication  systems,  a value  of  Pg  ■ 10  is  barely  acceptable;  any- 
thing larger  than  this  produces  undesirably  high  error  rates.  Consider 

the  "vacuum"  case,  S„  * 0 : Fig.  7.2  tells  us  t’ at  in  order  to  assure 

-5  1 

P £ 10  it  is  necessary  that  y >_  10.3  dB  . Still  using  this  figure, 

6 0 -5 

let  us  find  the  required  valte  of  y to  produce  P ■ 10  for  each 
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other  value  of  ; the  difference  between  this  yQ  and  10.3  dB  repre- 
sents the  additional  "vacuum"  signal-to-noise  ratio  necessary  in  order  to 
overcome  the  effects  of  fading  at  that  particular  level  of  scintillation. 
This  we  shall  call  the  link  degradation  loss  factor  L , which  is  plotted 
as  a function  of  S£  in  Fig.  7.3. 

Notice  the  appearance  of  a threshold  effect  in  this  figure,  repre- 
sented by  the  very  steep  rise  in  the  curve  beyond  a certain  point.  If 
we  arbitrarily  set  10  or  15  dB  as  the  maximum  acceptable  loss  factor  (the 
exact  value  is  not  very  critical),  it  is  clear  that  for  less  than 

about  0,25  the  system  can  operate,  but  for  S£  > 0.25  it  will  be  in  con- 
siderable difficulty. 


These  calculations  were  made,  of  course,  only  for  two  of  the  many 
different  kinds  of  communication  systems,  but  the  results  may  be  considered 
as  typical  of  the  effects  to  be  expected.  We  see  that  fading  due  to  scin- 
tillations causes  an  increase  in  the  probability  of  element  error  which 
may  reach  unacceptable  levels  if  the  link  signal-to-noise  ratio  is  too 
small  and/or  the  scintillation  index  is  too  large.  Because  the  natural 
environment  is  also  subject  to  slow  fading  whose  magnitudes  may,  with 
appropriately  small  probabilities,  cause  undesired  increases  in  P , 
satellite  communication  systems  are  usually  designed  to  have  a sufficiently 
large  (S/N)q  to  assure  satisfactory  performance  under  normal  fading  of 
given  probabilities  of  occurrence.  Figure  7.3  suggests  that  such  con- 
servatively designed  communication  systems  may  be  expected  to  survive 
artificial  scintillations  with  indices  up  to  approximately  S£  - 0.25  , 
but  it  is  unlikely  that  a nominal  system  will  have  been  designed  to  endure 
link  degradation  losses  in  excess  of  25  dB,  such  as  will  be  encountered 
when  > 0.3  . Accordingly,  we  can  say  that  scintillations  with  such 
high  indices  will  have  considerable  disruptive  effects  on  normal  communi- 
cation systems.  It  should  be  remarked,  however,  that  artificial  scintil- 
lations rarely  attain  these  high  values  of  S0  . 


APPENDIX 

SCINTILLATION  INDEXES  AND  STATISTICAL  DISTRIBUTIONS 


Several  quantitative  measures  for  the  scintillation  depth  of  signals 
from  radio  stars  and  satellites  are  in  common  use.  In  this  appendix  we 
describe  the  most  common  indexes  and  illustrate  some  of  the  interrelations 
among  them  for  some  models  of  the  scintillation  statistics. 

1/2 

In  measurements  of  satellite  signals  the  amplitude  R ■ (AA*) 

is  frequently  measured,  while  in  radio  star  measurements  the  power,  or  a 

2 

quantity  proportional  to  R , is  commonly  measured.  Briggs  and  Parkin 

have  defined  four  scintillation  indexes  based  on  the  mean  or  root-mean- 

2 

square  deviations  of  R or  R . These  ares 


|e(R2)  - [E(R)]2! 


E(R2) 


|e(R4)  - [E(R2)]2} 
s » J L 

4 E(R2) 

The  scintillation  index  measured  and  reported  by  Bowhill  and  others  at 
Aeronomy  Corporation  is  S£  (at  least  in  PRAIRIE  SMOKE  II-V) . Here 

is  equivalent  to  the  measure  S used  throughout  the  Briggs  and  Parkin 

2 2 
paper  and  is  equal  to  the  square  root  of  the  variance  of  R considered 

by  Me rcier. 


Preceding  page  blank 


95 


Given  the  probability  distribution  of  R , it  is  possible  to  find 
relations  among  S^,  S2»  , and  . However,  the  exact  distribution 

of  R is  not  generally  known.  Mercier  asserts  that  far  from  a random 
phase  screen,  R assumes  a Rice  distribution  with  density  function 

| exp[-(R2  + P2)/24']lo(RP/'t') 

in  Rice's  original  notation,  where  Iq  is  the  modified  Bessel  function 
of  zero  order.  The  parameter  ¥ is  proportional  to  the  amount  of  power 
in  the  signal  fluctuation;  P /2¥  is  a kind  of  "signal-to-noise"  ratio 
exactly  analogous  to  random  thermal  in-phase  and  quadrature  noise  added 
to  a signal  of  known  amplitude.  For  very  large  ¥ , the  fluctuations 
are  large  and  the  distribution  approaches  Rayleigh;  for  small  ¥ , the 
fluctuations  are  small  and  the  distribution  approaches  a shifted  Gaussian. 
Over  the  range  of  values  of  the  parameter  ¥ , the  scintillation  index 
as  a function  of  S^.  varies,  but  is  usually  about  0.5.  (The  ratio 
is  about  0.52  in  the  Rayleigh  limit.) 

29 

An  alternative  distribution  is  the  Nakagami  m-distribution  whose 
density  function  is  given  by 


2mm  R2m-1 

r(m)fim 


exp (-mR  /ft) 


2 

where  = E(R  ) . This  distribution  approximates  the  amplitude  of  the 

vector  sum  of  many  waves  with  random  amplitudes  and  phases.  For  large 

m the  distribution  approximates  the  Rice  distribution,  while  for  m ■ 1 

-1/2 

it  becomes  exactly  Rayleigh.  The  index  ■ m ; the  expression  for 
S2  is  moderately  complex,  but  Sj  - 0.5S^  as  before. 


Figure  A.l  shows  as  a function  of  S2  for  several  assumed 

distributions  of  R : Gaussian,  Rayleigh,  Rice,  and  Nakagami.  Except 


in  the  case  of  the  Rayleigh  distribution,  which  reduces  to  a single  point 
on  this  figure,  all  of  these  distributions  produce  curves  which  are  sub- 
stantially similar  to  ■ 2Sj  over  the  range  of  meaningful  values  of 
these  parameters. 

There  is  considerable  uncertainty  as  to  hich  distribution  best 

29 

fits  the  measured  data  from  natural  scintillation.  Nakagami  showed  good 

agreement  of  the  m-distribution  to  long-range  HF  (9.67  to  20.02  MHz)  com- 

30 

munication  channels.  Whitney,  et  al.  showed  fair  agreement  between 

the  Nakagami  distribution  and  the  distribution  of  fading  of  ATS-3  signals 

24 

at  136  MHz.  On  the  other  hand,  Rino  and  Fremouw  showed  that  the  Nakagami 
distribution  yields  only  a poor  approximation  to  ATS-3  fading  distribution 
measured  in  Peru.  They  indicate  good  agreement  to  a more  complicated  two- 
dimensional  Gaussian  (or  Hoyt)  distribution  for  which  the  in-phase  and 
quadrature  scattered  components  have  a joint  elliptical  one-sigma  contour. 
(The  Rice  distribution  is  a special  case  for  which  that  ellipse  is  a 
circle.)  This  distribution  can  yield  values  of  > 1. 

The  apparent  inconsistency  between  these  two  views  seems  to  stem 
from  different  sets  of  measurements.  Whitney  used  ATS-3  fading  data  mea- 
sured in  Massachusetts,  Greenland,  and  Peril.  The  best  fit  to  Nakagami 
statistics  was  observed  in  the  Massachusetts  data  for  which  the  fading 
was  less  severe  than  observed  in  Greenland  and  Peru.  Rino  and  Fremouw 
used  ATS-3  data  from  Peru  only,  with  stronger  fading.  The  Gaussian  model 
approaches  the  Rice  or  Nakagami  for  the  case  of  shallow  fading.  Conse- 
quently, the  apparent  inconsistency  seems  to  stem  from  the  different  range 
of  data.  We  should  note,  however,  from  the  results  of  Rino  and  Fremouw 
that  the  Nakagami  statistics  provides  a poor  fit  for  deep  fading. 

While  the  S-indices  and  the  m-index  can  be  computed  readily  from 
digitized  samples  of  the  signal  amplitude  or  power,  signal  levels  are 
frequently  recorded  only  on  paper  charts  so  that  these  exact  indexes  are 
difficult  to  determine-  Some  approximate  indexes  which  are  easy  to  deter- 
mine from  these  charts  are  in  common  usage.  One  form  is  to  compute 
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\ . (R2  “ Ri)/(R2  + Rp  * where  R2  and 

the  chart  in  some  way.  For  example,  R2 
f minimum  signal  levels  over  the  interval;  another  frequent  selection  is 

the  "third  peak  down,"  from  the  maximum  and  "third  peak  up"  from  the  mini- 

• 31 

mum.  Bischoff  and  Chytil  illustrate  the  relations  between  two  approxi- 
mate scintillation  indexes  ind  the  m-index. 

The  logarithm  of  the  signal  is  frequently  recorded,  especially  in 
conjunction  with  communication  satellites,  so  that  the  fading  can  be  ex- 
< pressed  in  dB.  The  dB  level  of  scintillations  is  commonly  the  difference 

in  dB  between  selected  peak  and  null  signal  levels  during  a recording 
interval.  With  the  Rice  or  Nakagami  distributions  for  signal  amplitude, 
almost  any  fading  depth  is  possible  with  some  small  probability.  However, 
one  can  relate  (or  m)  to  the  fading  depth  by  comparing  the  amplitudes 

at  two  specified  percentile  levels  of  the  distribution  function.  Figure 
A. 2 shows  the  range  in  dB  between  R^  and  Rq^,  and  R^  , and 

Rg9  9 and  Rq  ^ where  the  subscript  denotes  the  percentile  values  for 
the  Nakagami  distribution  of  amplitude. 


are  signal  levels  picked  from 
and  R^  can  be  the  maximum  and 
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